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ABSTRACT 

We report Hubhle Space Telescope/Cosmic Origins Spectrograph observations of the Lya emission 
and interstellar absorption lines in a sample of ten star-forming galaxies at 2 ^ 0.2. Selected on the 
basis of high equivalent width optical emission lines, the sample, dubbed “Green Peas,” make some of 
the best analogs for young galaxies in an early Universe. We detect Lya emission in all ten galaxies, 
and 9/10 show double-peaked line profiles suggestive of low H I column density. We measure Lya/Ha 
flux ratios of 0.5-5.6, implying that 5% to 60% of Lya photons escape the galaxies. These data confirm 
previous findings that low-ionization metal absorption (LIS) lines are weaker when Lya escape fraction 
and equivalent width are higher. However, contrary to previously favored interpretations of this trend, 
increased Lya output cannot be the result of a varying H I covering: the Lyman absorption lines (Ly/3 
and higher) show a covering fraction near unity for gas with Nh / > 10^® cm“^. Moreover, we 
detect no correlation between Lya escape and the outflow velocity of the LIS lines, suggesting that 
kinematic effects do not explain the range of Lya/Ha flux ratios in these galaxies. In contrast, we 
detect a strong anti-correlation between the Lya escape fraction and the velocity separation of the 
Lya emission peaks, driven primarily by the velocity of the blue peak. As this velocity separation is 
sensitive to H I column density, we conclude that Lya escape in these Green Peas is likely regulated 
by the H I column density rather than outflow velocity or H I covering fraction. 


1. INTRODUCTION 

The Lya emission line is a heavily-used diagnostic 
in studies of high-redshift galaxies and the intergalac- 
tic medium (IGM). At redshifts z > 2, this feature 
is shifted into readily observable optical wavelengths, 
leading many studies to rely on Lya for both discov¬ 
ery and spectroscopic redshift confirmation of galax¬ 
ies. Lya-based surveys are now extending to z ~ 6 
and beyond, where they are identifying the building 
blocks of present day galaxies and the sources that 
are most importa nt for metal enrichment and reioniza¬ 
tion of the IGM (iMartin et al.l 120081: lOuchi et al.l |2M^ 
Hu et al.l 120101: iKashikawa et al.l 12011 : Dressier et al l 


2011L 120151 IHenrv et al.ri2ni2l: IRhoads et al.1120121) " In 

addition, because Lya photons scatter resonantly in 
neutral hydrogen, they offer a wealth of information 
about this gas. Spectroscopic line profiles can con¬ 
strain gas content and kinematics ( Verhamme et all 


12014lErb et al.l[20l4 IMartin et al.ll2015ll . and Lya imag¬ 

ing can illuminate the gaseous halos a r ound galax¬ 
ies (Zhen g et al.l 120111: ISteidel et all 120111 iHaves et al.l 
I20H: : iMomose et al.l I2014IL In the epoch of reion za- 
tion, weaker or less frequent Lya emission may indi- 
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cate a significant neutral hydrogen frac t ion in the IGM 
(lota et al.l 120101 iPentericci et al.l 120111 : iSchenker et al.l 
120121: iTreu eTall 120121 I20lli , an incre ase in the es¬ 
cape fraction of ionizing photons, or both (iDiikstra et al.l 

[Ml . 

Despite the fact that a wide range of studies rely on 
Lya emission, we do not understand how to interpret the 
feature. The Lya luminosity, equivalent width, ^Lya, 
and Lya/Ha flux ratio of galaxies are seen to vary widely, 
and are difficult to predict from other galaxy proper- 
ties like star formation rate (SFR) and dus t extinction 
(|Giavalisco et al.l 119961: IScarlata et al.l I2009fl . Qualita¬ 
tively, it is believed that the Lya emission line is mod¬ 
ified by resonant scattering in the interstellar medium 
(ISM) and circumgalactic medium (GGM). This process 
increases the optical depth of Lya photons, making them 
more susceptible to dust absorption and sometimes turns 
the emission profile into absorption. Alternatively, in the 
absence of dust, the Lya photons may simply diffuse to 
large radii and be misse d by poor sensitivity to low sur¬ 
face b rightness emission (ISteidel et al.ll20lil : IHaves et al.l 
IMl . Besides dust and geometry of the HI gas, the 
kinematics of galaxy outflows may also play a role, since 
Lya photons may escape more easily when scattered by 
neutral hydrogen that is out of resonanc e with the ISM 
(|Kunth et al.lli99§ : iWofford et al.l 120131) . Nevertheless, 
our poor quantitative understanding of Lya output is 
cause for concern when studies try to draw robust con¬ 
clusions about galaxies and the IGM from this feature. 

Efforts to understand Lya escape from high red- 
shift galaxies have generally suffered from the limita- 
tions associated with studying faint, distant sources (e.g. 
Shaplev et al.l 120031: lErb et ahi I2010L 12014 iJones etahi 
2012L 1201311 Fortunately, many of these difficulties 
can be overcome when studying nearby galaxies. If 
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TABLE 1 

Green Pea Sample & Observations 


ID 

RA 

(J2000) 

DEC 

(J2000) 

Z 

E(B-V)mw 

(mag) 

Rest Wavelength Coverage 

(A) 

G130M Exposure 

(s) 

G160M Exposure 

(s) 

0303-0759 

03 03 21.41 

-07 59 23.2 

0.164887 

0.0877 

975 - 1515 

2190 

3829 

1244+0216 

12 44 23.37 

02 15 40.4 

0.239420 

0.0213 

945 - 1430 

2042 

6507 

1054+5238 

10 53 30.80 

52 37 52.9 

0.252645 

0.0132 

910 - 1425 

824 

2736 

1137+3524 

11 37 22.14 

35 24 26.7 

0.194396 

0.0161 

965 - 1505 

1264 

2340 

0911+1831 

09 11 13.34 

18 31 08.2 

0.262236 

0.0248 

900 - 1435 

2074 

6530 

0926+4427 

09 26 00.44 

44 27 36.5 

0.180698 

0.0165 

970 - 1505 

5640 

6180 

1424+4217 

14 24 05.72 

42 16 46.3 

0.184800 

0.0094 

965 - 1220 

1209 

Oa 

1133+6514 

11 33 03.80 

65 13 41.4 

0.241400 

0.0097 

945 - 1430 

1232 

4589 

1249+1234 

12 48 34.63 

12 34 02.9 

0.263403 

0.0252 

900 - 1425 

1644 

6372 

1219+1526 

12 19 03.98 

15 26 08.5 

0.195614 

0.0239 

965 - 1505 

716 

2304 


® Failed Observation 


low-redshift galaxies can be found where the physi¬ 
cal conditions are similar to high-redshift star-forming 
galaxies, it is possible to learn more about the astro¬ 
physics that regulates Lya output. Despite this re¬ 
alization, most investigations to date have struggled 
to identify local galaxies that are similar to the low- 
mass, low-metallicity, high specific SFR galaxies that 
dominate the cosmic star-formation at early times (e.g. 
lAlavi et all 12014 IStark et al.l 120141 . Instead, present 
samples of nearby galaxies with Lya observations are 
on average, more massive, dustier, metal-rich, and 
have lowe r specific SFRs than th e bulk of high-r e dshift 


Heckman et al. 2011: Cowie et al.l 12011: Wofford et al. 

201 IHaves et aHI201.81120141: IMartin et al.ll201.'in. As a 

possible consequence of 
in low redshift samples- 
order of magnitude low 
Lva-selected galaxies of 

this selection, Lya luminosities 
when measured- are around an 
3r than the luminosities of most 

iserved at z > 3 (IGronwall et al. 

2007 

: lOuchi et al.l 12010 

: IHu et al.l 120101: iDressler et al. 

2011 

, 120151: iHenrv et al. 

120121) 


The recent discovery of extremely high equivalent 
width emission line galaxies (the “Green Peas”) within 
the Sloan Digital Sky Survey (SDSS) offers a new avenue 
to investigate Lya emis sion from high-redshift analogs 
(|Cardamone et ahl 120091) . With masses reaching below 
10® Mq, metallicities lower tha n 12 -|- log(0/H) ^ 8.2 
(from direct Te measurements; lizotov et al.l 1201 ill , and 
Ha equivalent widths {WHa) exceeding hundreds of A, 
these local galaxies may be more representative of high 
redshift Lyman Alpha Emitters (LAEs; 20 A) 

and reionization epoch galaxies than other local samples. 
As a comparison, inferred emission line contamination 
to Spitzer/IRAC photometry provides tentative evidence 
for similar W Ha i n star-forming ga l axies at z 4 — 7 
(iShim et al.l 120111: IStark et all 120131 : iLabbe et al.l 120131 : 
iSmit et al.l 1201411 . Ultimately, the James Webb Space 
Telescope will clarify which segment of the high-redshift 
populations are most analogous to the Green Peas. At 
present, however, quantifying the Lya output from these 
nearby objects is a critical benchmark for future compar¬ 
isons. 

In this paper we present a UV spectroscopic study 
of Green Peas using the Cosmic Origins Spectrograph 
(COS) on the Hubble Space Telescope (HST). For a sam¬ 
ple of ten galaxies at 0.16 < z < 0.26, we measure the 
Lya emission that escapes the galaxies from within a 
few kpc (corresponding to the COS aperture). At the 


same time, by observing absorption lines in both hydro¬ 
gen and metals, these data allow us to quantify the role 
of galactic outflows, ISM, and CGM gas. This paper is 
organized as follows: in §2 we describe our sample se¬ 
lection and the galaxy properties derived from GALEX 
and SDSS; in §3 we describe the COS observations and 
resultant data and in §4 we present the Lya line profiles 
and measurements; §5 explores how the stellar popula¬ 
tion properties and dust of the Green Peas compare to 
other nearby samples that have been observed in Lya. 
In §6 and §7 we show how the strength of the Lya line 
varies with interstellar absorption equivalent width and 
kinematics, and discuss the velocity structure of the Lya 
line. Finally, §8 offers a comprehensive interpretation of 
the data and our conclusions are s ummarized in jjQ . 

Throughout this paper we use a IChabrieil (|2003f) ini¬ 
tial mass function (IMF), and a ACDM cosmology with 
Dm = 0.3, Ha = 0.7, and Hq = 70 km s“^ Mpc“^. We 
adopt the equivalent width sign convention used in liter¬ 
ature focused on high-redshift galaxies: positive equiv¬ 
alent widths indicate emission while negative values are 
used to signify absorption. All COS, GALEX and SDSS 
data are corrected for Milky Way extinction using at- 
ten uation measured bylSchlaflv fc Finkbeineil (1201111 and 
the iFitzpatrickl (1199911 extinction law. The foreground 
extinction values were obtained from the NASA Extra- 
galactic Databas^H, and are listed in Table [TJ 

2. SAMPLE SELECTION AND PROPERTIES 

The present sample of Green Peas were drawn from 
the catalog presented bv iCardamone et al.l (1200911 . To 
systematically id entify objects origi nally discovered by 
the Galaxy Zoo (|Lintott et al.lf2008ll . Cardamone et al. 
defined a color selection for objects with strong [O III] 
emission in the SDSS r— band. Following this selection, 
objects with low S/N SDSS spectra were removed and 
AGN were rejected on the basis of broad emission lines 
or t heir optical emission line ratios (Ha/[N II] vs. [O III]/ 
Hd: lBaldwin et anil981jl . The remaining sample consists 
of 80 galaxies with rest frame [O III] A5007 equivalent 
widths of hundreds to over 1000 A. 

We selected Green Peas that were bright enough in the 
FUV to enable continuum dete ction with COS. From th e 
80 star-forming Green Peas in ICardamone et all (|2009ll , 
we considered objects with GALEX photometry (GR6) 
and niFuv < 20 AB. We also required z < 0.27, so 

^ http://ned.ipac.caltech.edu/ 
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TABLE 2 

Green Pea Properties from GALEX+SDSS 


ID 

(1) 

Who. 

(A) 

(2) 

Ha/H/3 

(3) 

E{B - V)gas 
(mag) 

(4) 

logfDffa/ergs U 

(5) 

SFR 

(M© yr-i ) 
(6) 

Log M/M0 

(7) 

12 -h log(0/H) 

(8) 

Mpuv 

(mag) 

(9) 

d 

(10) 

Ruv 

(kpc) 

(11) 

0303-0759 

670 

2.78 

0.00 

42.24 

7.6 

8.89 

7.86 

-20.35 

-2.23 

0.8 

1244-1-0216 

840 

3.10 

0.07 

42.78 

26.2 

9.39 

8.17 

-20.32 

-1.70 

2.6 

1054-1-5238 

400 

3.15 

0.08 

42.71 

22.4 

9.51 

8.10 

-21.31 

-1.94 

1.3 

1137-1-3524 

580 

3.08 

0.06 

42.58 

16.8 

9.30 

8.16 

-20.56 

-1.78 

1.8 

0911-1-1831 

420 

3.50 

0.17 

42.68 

21.1 

9.49 

8.00 

-20.56 

-1.82 

1.1 

0926-1-4427 

610 

3.20 

0.10 

42.49 

13.6 

8.52 

8.01 

-20.58 

-1.98 

1.0 

1424-1-4217 

1100 

3.01 

0.04 

42.57 

16.5 

8.08 

8.04 

-20.40 

-1.90 

1.0 

1133-1-6514 

300 

2.90 

0.01 

42.02 

4.6 

9.04 

7.97 

-20.40 

-1.93 

1.9 

1249-1-1234 

670 

3.09 

0.07 

42.50 

13.8 

8.79 

8.11 

-20.25 

-1.82 

1.8 

1219-1-1526 

1270 

2.87 

0.00 

42.43 

11.9 

8.09 

7.89 

-19.94 

-1.65 

0.7 


Note. — Quantities are derived from the SDSS and GALEX, or taken from fT^tov et ^ l)2QlUh Measurements errors are dominated 
by systematics in these high S/N data, so statistical errors are not given. Emission line measurements include both broad and narrow 
components, as described in J2l Column Descriptions: (1) Object ID (2) Rest-frame Ha equivalent width, in A, measured from SDSS 
spectrum. (3) Flux ratio, F(Ha) / F(H/3) measured from SDSS spectrum (4) Nebular gas extinction inferred from the Hq/H/ 3 flux ratio, 
assuming the intrin sic ratio of 2.86 a s described in the text. (5) Ha luminosity measured fr om the SDSS sp ectrum (6) SFR estimated 
from Lffn u s ing th e IKennicuttj d 19981) calibration, di vided by 1 .8 to convert from a Salpeter to lChabrieil d2003i) IMF (7) Stell ar mass from 
I Izotov et al.l (|2Q1 II) . also converted from a Salpeter to lChabrierl (|2003l ) IMF (8) Direct, Te([0 Hl])-based metallicities taken from lizotov et aTI 
J2011I) (9) FUV absolute magnitude (AB), calculated from GALEX photometry, interpolated to provide an estimated rest-frame 1500A 
luminosity (10) UV slope, defined as Fx oc A^, calculated from GALEX data. (11) NUV Petrosian radius (see text in la, calculated from 
the COS acquisition image. 


that the FUV spectra would cover the Si IV AA 1393, 
1403 lines. The resulting sample of ten galaxies is listed 
in Table [H We verified that these galaxies have stel¬ 
lar masses, SFRs, metallicities, and [O III] equivalent 
widths consistent with the parent sample of Green Peas 
(|Cardamone et ahl 120091 lizotov et ^ 1201111 . However, 
because of the FUV magnitude limit, the average UV 
luminosity in the present sample is around 1.6 times 
hi gher than the average (G ALEX detected) Green Peas 
in iCardamone et ^ (| 20091) . 

In Table d] we compile properties of the Green Peas, 
drawn from the literature or derived from GALEX and 
SDSS. We describe these quantities below. In most cases, 
the high S/N of these data imply insignificant statistical 
errors on the measurements. Therefore we do not list 
these uncertainties. 

First, we chose to re-measure the emission lines in 
the SDSS spectra rather than use catalog measurements. 
One of the Green Peas in our sample (0303-0759) con¬ 
tains an unphysical Balmer decrement (an Ha/H/3 flux 
ratio of 2.25) in the MPA-JHU DRi(0 catalogs, which is 
the result of an incorrect Ha flux. This failure may be at¬ 
tributable to difficulty in fitting the nearly blended Ha-|- 
[N II] lines with single Gaussians, when the lines clearly 
show evidence of broad wings extendin g h undreds of km 
s~^ ( also noted bv lAmorfn et al.l[2012l and iJaskot fc OevI 
120141) . Therefore, for each of the Green Peas in this pa¬ 
per, we re-fit the emission line spectra with two kinematic 
components. We defined a kinematic model with a nar¬ 
row component centered at the systemic redshift (u = 0), 
and a broad component which is allowed to be offset from 
u = 0. For 25 of the strongest lines in the spectrum (Ha 
through H9, [S II], (N II], (O III], [O II], (Ne III] and 7 
He I lines), the amplitudes of the broad and narrow com¬ 
ponents were allowed to vary, but their velocity widths 
and centroids were required to be the same from line 
to line. Hence, the resultant model has 54 free param- 

® http://www.mpa-garching.mpg.de/SDSS/DR7/ 


eters (redshift, Doppler shift of the broad component, 
velocity widths of both the broad and narrow compo¬ 
nent, and two Gaussian amplitudes for each line). This 
model was fit to the continuum subtracted spectrum pro¬ 
vided by the SDSS, thereby measuring redshifts (Table 
[T]), emission line fluxes (for the combined broad and nar¬ 
row kinematic components), and the kinematics of the 
ionized gas. Equivalent widths are evaluated by taking 
the ratio of the emission line flux to the median contin¬ 
uum flux density in a 100 A window centered on each 
line. The kinematic measures derived from our fits are 
similar for all of the Green Peas, showing a narrow com¬ 
ponent at the systemic redshift with a FWHM around 
200 km s“^, and a broad component, also centered near 
u ^ 0, with a FWHM around 500 km s“^. For the iso¬ 
lated H/3 line, the broad component contributes 15-35% 
of the total line flux in these galaxies. 

Table [2] lists the Ha equivalent width and the Ha to 
H/3 flux ratio. From this latter quantity, we derived 
dust extinction assuming a lGalzetti et al.l ()2000f l extinc¬ 
tion curve and an intrinsic ratio of F(Ha)/F(H/3) = 2.86. 
This intrinsic ratio is appropriate for electron tempera¬ 
ture Te = 10^ K and density = 10^ cm“^; increasing 
Te two-fold decreases the intrinsic Balmer decrement to 
2.75 and increases the extinction, E{B — V)gas, by ap¬ 
proximately 0.03 magnitudes. Dust corrected Ha lumi¬ 
nosities are given, and converted to SFR usin g the cali¬ 
brati on given in Kennicuttt (1998; adjusted to a lGhabri^ 
120031 initial mass function). 

The stella r mass es of the Green Peas are taken from 
lizotov et al.1 ()2011h . These authors note that, for high 
equivalent with emission lines, the corresponding bound- 
free continuum emission can be comparable to the stel¬ 
lar continuum. Therefore previous estimates of stellar 
mass in these galaxies may be too high; to correct for 
this mis-estimation, Izotov et al. derived stellar masses 
from fits to the SDSS spectra, using models that include 
nebular continuum. For consistency, we have decreased 
these masses by a factor of 1.8 to convert them from a 
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Salpeter to iChabrieil (|2003D initial mass function. Oxy - 
gen abundances are also taken from lizotov et al] (1201 If) , 
where they were derived from [O III] electron tempera¬ 
ture ([O III] A4363). All 10 of the Green Peas in this 
paper show detections in this line, confirming their low 
metallicities. 

The ultraviolet continuum luminosity and slope were 
derived from GALEX photometry (FUV: A ^ 1340 — 
1790 A, NUV: A - 1770 - 2830 A). The FUV fluxes 
were corrected for the contribution from Lya emission of 
the Green Peas (described in ]j4]), and both bands were 
corrected for Galactic foreground extinction (as noted in 
§1). From the corrected FUV - NUV color, we calculated 
the UV power law slope, f3, {Fx oc A^]0- The absolute 
FUV magnitude was calculated by interpolating between 
the FUV and NUV bands to obtain the luminosity, L^, 
at 1500A. 

In the final column of Table [2] we list Ruv, the size 
of the galaxy measured from the NUV acquisition im¬ 
ages. For the purposes of compariso n in j|5l we choose a 
Petrosian radius, rp 2 o as adopted bv iHaves et al.l (I2014D . 
This radius defines the circular isophote where the local 
surface brightness is 20% of the internal surface bright¬ 
ness. 

3. DATA 

3.1. Observations 

The observations in this paper are part of two GOS 
programs. One of the galaxies that met our selection 
(0926-1-44 27) was also classified as a Lyman Break Ana¬ 
log (LBA: [Heckman et al.ll20lil) . and was previously ob¬ 
served with COS (GO 11727; PI T. Heckman). The re¬ 
maining 9 Green Peas were observed as part of GO 12928 
(PI A. Henry). 

For all ten of the galaxies, the target acquisition was 
accomplished with NUV imaging, configured with the 
Primary Science Aperture (PSA) and Mirror A. As de¬ 
tailed in the COS Instrument Handbook, an initial im¬ 
age is obtained and analyzed to find the peak NUV flux. 
Next, HST is moved to place the peak flux in the cen¬ 
ter of the COS aperture, and a second image is taken 
to verify the shift. The target acquisition images of the 
Green Peas are shown in Figure |T| These data indicate 
relatively compact emission, implying that most of the 
UV continuum emission from these galaxies falls within 
the central 1" diameter of the COS aperture (shown as 
circles) where vignetting is minimal. (Likewise, if opti¬ 
cal line emission closely follows the UV continuum, we 
expect negligible aperture losses in the 3" SDSS fiber 
spectroscopy.) 

The spectra were obtained in the FUV, using both the 
G130M and G160M configurations (in 9/10); with both 
gratings the rest-frame wavelength coverage spans from 
approximately 950 — 1450 A, or somewhat shorter/longer 
depending on the redshifts. The rest-frame wavelengths 
covered for each galaxy are listed in Table lU along with 

® The UV slope has been s hown to differ syste matically depend¬ 
ing on the adopted bandpass (IHowell et aLll^lOll . Therefore, these 
UV slopes are not strictly comparable to 0 derived from the In¬ 
ternational Ultraviolet Explorer (lUE; e.g ICalzetti et aTlll994ll . or 
GALEX observations of galaxies at z ~ 0. We take the present 
measurement as a reasonable approximation; our conclusions do 
not rely on /?. 


the exposure time for each grating. For 1424-1-4217, the 
second orbit containing the G160M observations failed; 
it did not qualify for a repeat observation because the 
program was more than 90% complete. In this case, the 
observations cover the Green Pea Lya emission only for 
velocities v < 750 km s“^. As recommended by the 
COS Instrument Handbook, all spectra were taken at 
four FP-POS settings in each grating in order to miti¬ 
gate the effects of fixed pattern noise. In addition, the 
central wavelengths for each grating were chosen to avoid 
placing strong absorption lines between the A and B seg¬ 
ments. For one galaxy (1054-1-5328) this choice was not 
possible in G130M, so two central wavelengths were used. 
Since CALCOS does not combine data taken at different 
central wavelengths, we combined them by simply tak¬ 
ing the mean (except in the gaps, where we adopt the 
values from the spectrum that has coverage). The errors 
were derived from taking the Poisson noise on the total 
counts, as we will describe below in (13.31 

All data from GO 12928 were taken at the COS life¬ 
time position two; they were downloaded from MAST 
2013 August 23 implying that they were processed with 
CALCOS version 2.19.7. The spectra from GO 11727 
were taken in lifetime position one, and were downloaded 
on 2012 November 2 (CALCOS version 2.18.5). 

The full spectra for all ten Green Peas are shown in 
the Appendix. 

3.2. Estimating Spectral Resolution 

Because the Green Peas are spatially resolved, the 
spectral resolution that we achieve is reduced from the 
15 - 20 km s“^ that is expected for point sources. Nev¬ 
ertheless, Figure [T] shows that the galaxies are compact; 
they do not fill the aperture, so the resolution should 
be better than the 200 km s~^ that is rep orted for filled 
aperture observations ([France et al.ll200^ . 

In order to quantify the spectral resolution for the 
present sample, we create a model line spread func¬ 
tion (LSF) from the NUV acquisition images and the 
point source LSPF^. For each galaxy, we create a one¬ 
dimensional profile in the dispersion direction, by sum¬ 
ming the pixels in the cross dispersion direction of the 
acquisition image. This profile is then convolved with the 
LSF to estimate how it is broadened by spatially resolved 
emission. 

The resultant model LSFs in both G130M and G160M 
have typical FWHM ^12 pixels (or 18 for the somewhat 
larger 1244-1-0216) at 1450A. This scale is larger than 
the 6-8 pixel FWHM of the point-source LSF, although 
both the point-source and model galaxy LSFs have wings 
extending beyond ±20 pixels. For a dispersion of 9.97 
(12.23) mA pixel-1 in G130M (G160M), the FWHM of 
the LSF corresponds to 25 (30) km s”! for typical Green 
Peas and 37 (46) km s”! for 1244±0216. As will be¬ 
come evident in the later sections, the absorption lines 
are much broader and clearly spectrally resolved. The 
Lya emission may have somewhat poorer spectral reso¬ 
lution if it is more extended than the stellar component. 

3.3. Binning and Noise 

10 


http://www.stsci.edu/hst/cos/performance/sp6ctral_resolution/ 
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Fig. 1. — COS/NUV acquisition images are shown for the present sample of Green Peas. Cutouts are 1.6^' on a side, shown oriented 
with north up and east towards the left. The circle, with a 1" diameter, approximates the unvignetted portion of the COS aperture, while 
the full 2.5^^ diameter aperture extends beyond the edge of the images. The arrows indicate the dispersion direction in the spectra. The 
images are displayed in order of increasing Lyo; escape fraction (left to right, top to bottom; given in Tableland described in S- 


We chose to bin the spectra to gain signal-to-noise 
and mitigate the p otential effects of correlated noise 
(|Keenev et alJ 1^121 . As we showed above, the spectral 
resolution for our galaxies is degraded relative to our ex¬ 
pectation for point sources. Binning 20 native COS pix¬ 
els into one does not discard information for the galaxies 
in our sample. Therefore, we used the IDL rebin func¬ 
tion to calculate the mean of 20 adjacent fluxes, f\ and 
wavelengths, A. 

Inspection of the error arrays for the unbinned spectra 
showed that CALCOS overestimated the statistical un¬ 
certainties in the nine Green Peas observed in GO 12928, 
which have lower S/N than the spectrum of 0926-l-442'<F^. 
We measured the noise in line-free regions of these spec¬ 
tra, and found that it was around three times smaller 
than implied by the error vector. Therefore, we re¬ 
calculated the error spectrum by taking the Poisson noise 
on th e observed cou nts in each 20 pixel bin, accord¬ 
ing to iGehrelsI l|1986ll . The total counts per pixel were 
taken from the GCDUNTS array in the one-dimensional 
COS spectrum, which includes both source and back¬ 
ground counts. To convert the error on the counts to 
an error on the flux, we divide by the exposure time 
and the sensitivity function. The sensitivity function 
was simply estimated from the ratio of the calibrated 
and uncalibrated one-dimensional spectra (the FLUX and 
NET arrays provided by CALCOS). We verified that the 
newly derived error vector is a sensible representation of 
the noise in line-free portions of the continuum. More¬ 
over, this method applied to the higher S/N spectrum of 
0926-1-4427 showed good agreement with the CALCOS 
error spectrum. For the sake of uniformity with the rest 
of the sample, we adopt the re-calculated error spectrum 
for this object. 

Significant amounts of non-Poissonia n (correlated) 
noise have been noted in COS spectra (|Keenev et al.l 

This mis-estimation is a known problem for low S/N data 
(Oliveira 2014, private communication). 


[Ml), so we also test for this effect in our data. For 
this analysis, we repeated the spectral binning described 
above, creating spectra with V = 3, 5, 10, 20, 40, 100, 
and 200 COS pixels averaged. We also include an un¬ 
binned spectrum to probe N = 1. In each spectrum, 
we subtracted a local continuum fit and measured the 
noise, (Jn, in regions that were free from strong absorp¬ 
tion lines. In the case of the Poisson limit, we expect the 
noise per bin to scale with the noise per single COS pixel 
(Tat = uiN^, where (3 = —0.5. In reality, we measure 
/3 = —0.43 to —0.48, close to the Poisson limit. This 
contribution from non-Poissonian noise is smaller than 
the contribution measured by Keeney et al; it does not 
impact our conclusions so we do not correct for it. 

3.4. Velocity Precision 

Since we are interested in kinematic features probed 
by our COS spectra, it is important to assess the wave¬ 
length solution and velocity precision of our data. We 
make three tests. First, we check that the geocoronal 
Lya 1215.67 and O I A1302.17 lines lie at the correct 
observed wavelengths in the G130M spectra. To make 
this comparison, we first remove the heliocentric velocity 
correction implemented by CALCOS. This step ensures 
that the geocoronal emission is in the rest-frame of the 
Earth. After taking this correction into consideration, 
the geocoronal emission features fall within -5 to -|-17 
km s“^ of their expected velocities. 

The second test that we make is a comparison of the 
Milky Way ISM absorption line velocities. Because these 
lines may be Doppler shifted, they do not test the zero- 
point accuracy of the wavelength scale. Rather, their 
consistency gives a measure of the precision of our ab¬ 
sorption line velocities. The features that we use for 
this comparison are the Si II AAII90.4, 1193.3, AI260.4, 
and AI526.7, as well as C II A1334.5, and Al II AI670.8. 
Although higher ionization states are sometimes de¬ 
tected, we do not include these lines becuase they can 
have different kinematics than the low-ionization lines 
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Fig. 2.— The 10 Green Pea galaxies observed with COS all show 
Lya in emission, even though they were not selected to exhibit this 
feature. Broad wings on the lines clearly extend to several hundred 
km s“^, and 9/10 show double-peaked emission. The profile of 
1424-1-4217 is truncated at the red-cutoff of the G130M spectrum; 
for this object, the G160M observation failed. The spectra are 
displayed in order of increasing (left to right, top to bottom). 


(when observed in other galaxies, e. g. IGrimes et al.ll200^ 
iSteidel et al.ll201(il) . From Gaussian fits to the Milky Way 
lines, we find velocities that are consistent within 20-40 
km s“^. 

Finally, we confirm that the systemic redshift of the 
stars in the Green Peas are consistent with the H II re¬ 
gions. Redshifts measured from the photospheric C III 
1175.5 line are compared to the redshifts that we mea¬ 
sured from the SDSS spectra (Table[T|). In 8/10 galaxies, 
these redshifts agree to better than 41 km s“^; for two 
others, the test is not possible: the G III line is unde¬ 
tected in 1219-1-1526 and it is contaminated by Milky 
Way absorption in 1137-1-3524. Nevertheless, when mea¬ 
surable, the velocity offset between the H II regions and 
stars are consistent within the uncertainties. 

In summary, these tests show that we can measure 
velocities in the present GOS FUV spectra to better 
than 40 km s“^. This level of precision is confirmed for 
COS observations o f Ultra Luminous Infrared Galaxies 
(jMartin et al.l[20T^ . 

4. Lyo EMISSION FROM GREEN PEA GALAXIES 

4.1. The prevalence of strong, double peaked Lya 
emission 


Figure [5] shows the Lya spectra of the Green Peas 
in the present sample, displayed in order of increasing 
Lya escape fraction (estimated below). Their appear¬ 
ance is remarkable, for a few reasons. First, all ten 
of the galaxies show Lya in emission. This result is 
not trivial; typical star-forming galaxies at all redshifts 
show a range of Lya strength, ranging from pure (even 
damped) absorpt ion to P-cygni emiss i on plus absorption , 
to pure emission (iWofford et alJl2013HPettini et al.ll2002l: 

iShaolev et al.ll2003D . On the other hand, the Green Peas 
are not typical of nearby star-forming galaxies, so the 
Lya emission seen in all ten galaxies is not entirely sur¬ 
prising. Nearby, high equivalent width optical emission 
lines, like those in the Green Peas, have been suggested 
to favor Lya emission (|Cowie et al.ll201ll) . Furthermore, 
in more distant galaxies at 4 < z < 6, the increasing in- 
cidence and strength of Lya emission (IStark et al.ll20Hl) 
may go hand in hand with high equivalent width opti¬ 
cal emission lines (i nferred from contamination to broad¬ 
band photometry ; iShim et all 120111 : IStark et all 120131 : 
iLabbe et all 1201^ . The detection of Lya emission in 
all ten Green Peas is consistent with these claims. 

The second noteworthy feature about the spectra in 
Figure [2] is the prevalence of the double-peaked line 
shape. Of the Green Peas that we observed with 
COS, 9/10 share the spectral-morphology with both red- 
shifted and blue-shifted emission peaks. Again, this line 
sha pe is not typica l of nor mal nearby star-forming galax¬ 
ies. IWofford et al.l (|2013D find only one doubled-peaked 
emitter in their sam ple of twenty Ha se lected objects 
around z ^ 0.03. And iMartin et alJ (I2015D . in a study of 
eight nearby ULIRGs, show that Lya can exhibit com¬ 
plicated kinematic profiles with broad blueshifted emis¬ 
sion in more dusty envi ronments. In contr ast, at higher 
redshifts (z 2 — 3), iKulas et al.l (j2012[) report that 
30% of UV-continuum selected galaxies with Lya emis¬ 
sion show multiple-peaked profiles. Comparison to the 
Green Peas, however, reveals the importance of spec¬ 
tral resolution in this measurement . At th e 200-500 km 
s“^ resolution used bv IKulas et al.l (|2012D . many of the 
Green Pea spectra would be observed as single-peaked 
lines. Indeed, spectroscopy with 120 km s“^ resolution 
sh owed double-peaked line prohles in 3/3 LAEs targeted 
bv iChonis et al.l (|2013l) . Nevertheless, the Lya line pro¬ 
files of the Green Peas are notably different than the 
other low-redshif t samples, where ^ Hr, is lower and dust 
content is higher ( IWofford ed al.l[MT3 : lMartin et al.ll2015l : 

iRivera-Thorsen et al.11201,^ As indicated by Lya ra¬ 

diative transfer models, the blue peak should appear 
when the H I column density is low, plausibly due to 
anisotropies in the gas distribution (iVe rhamm e et al.l 
12014 [Behrens et al.ll20l4 iZheng fc Wallacell 20 141 ). 

In order to better understand how the unusual con¬ 
ditions in the Green Peas are influencing their Lya es¬ 
cape, we next provide quantitative measurements of the 
Lya lines and show how these compare to nearby sam¬ 
ples. Then, in the remainder of the paper we will use the 
UV absorption lines to explore how the conditions in the 
ISM and CGM affect the Lya emission. Radiative trans¬ 
fer models of the Lya spectral line profiles will follow in 
Orlitova et al. (in prep). 

4.2. Lya Emission Line Measurements 
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TABLE 3 

Lyo measurements from cos Spectroscopy of Green Pea Galaxies 


ID 

^Lya 

(10-14 gj.g g-1 cm“^) 

^Lya 

(10'^^ erg s“^) 

Lya/Ha 

rLyoi 

J esc 

^^Lya 

U/red 

Lya 

(A) 

•f^blue 

Lya 

^veak 

%ed 

,.peak 

^blue 

(km 

^,red 

^max 

s-l) 

^.blue 

^max 

0303-0759 

1.1 ±0.2 

0.8 

0.5 

0.05 

9±2 

9±2 

-0.1 ± 0.3 

170 

-290 

900 

-400 

1244+0216 

2.0 ±0.1 

3.4 

0.7 

0.07 

48 ± 10 

36 ±7 

12 ± 2 

250 

-280 

1100 

-1000 

1054+5238 

1.7 ±0.2 

3.1 

0.8 

0.07 

12 ±3 

12 ±2 

1±1 

160 

-250 

700 

-700 

1137+3524 

3.8 ± 0.2 

4.0 

1.2 

0.12 

35 ±7 

33 ±7 

3±1 

150 

-400 

1100 

-800 

0911+1831 

3.3 ±0.1 

6.8 

2.3 

0.16 

59 ±12 

48 ± 10 

8±2 

90 

-280 

1100 

-800 

0926+4427 

6.0 ±0.3 

5.4 

2.3 

0.20 

40 ±8 

36 ±7 

5±1 

160 

-250 

1200 

-800 

1424+4217 

8.5 ±0.2 

8.0 

2.4 

0.25 

95 ± 19 

57 ±12 

35 ± 7 

230 

-150 


-850 

1133+6514 

2.1 ±0.1 

3.6 

3.5 

0.40 

36 ±7 

25 ±5 

11 ±2 

230 

-100 

1000 

-600 

1249+1234 

5.4 ±0.1 

11.3 

4.4 

0.41 

98 ±20 

74 ±15 

16 ±3 

70 


1300 

-700 

1219+1526 

13.7 ±0.2 

14.7 

5.5 

0.62 

164 ± 33 

118 ±24 

41 ±8 

140 

-100 

1300 

-950 


Note. — Lya measurements from the COS spectra of the Green Peas give the fraction of Lyo that escapes from within the few kpc probed 
by the spectroscopic aperture. The Lyo fluxes and equivalent widths are calculated by directly integrating the line profile to the velocity where 
it meets the continuum, and and calculated for v > 0 and v < 0, respectively. The fluxes and equivalent widths assume 

a conservative 20% error on the continuum flux. Additionally, Lyo kinematics give the velocities marking the red and blue peaks, as well as the 
maximal blue and red velocities where the emission reaches the continuum. The former velocities are good to better than 50 km s“^ in all cases 
except for the red peak of 0926+4427, which appears impacted by foreground absorption. The maximal velocities a less certain, with typical 
errors around 200 km s“^. 


The Lya emission line measurements are presented in 
TablelH The line flux is measured by directly integrating 
the emission line spectra out to the velocity where the 
continuum is met. We adopt a conservative 20% uncer¬ 
tainty on the continuum flux near Lya, since the broad N 
V 1240 A stellar wind feature makes continuum normal¬ 
ization difficult. This error makes little difference for the 
Lya flux uncertainties, but it dominates the error bud¬ 
get for equivalent width measurements. For 1424-1-4217, 
where the G160M spectrum was unavailable, we use the 
G130M observations even though they only cover Lya 
for V < 750 km s“^. Since the other galaxies show only a 
small contribution at these velocities, we do not correct 
Lya measurement for missed flux. 

Although the compact sizes of the Green Peas sug¬ 
gests minimal aperture losses in the continuum spectra, 
the Lya emission may be more extended. At the red- 
shifts of our sample, the unvignetted portion of the COS 
aperture shown in Figure[T](l" diameter) corresponds to 
2.8 - 4.0 kpc; the full COS aperture is 2.5 times larger. 
As a comparison, in the nearby Lyman Alpha Refer¬ 
ence Sample (LARS), which is comprised of 14 galax¬ 
ies with HST imaging of Lya, all but two objects show 
emission extending to at least 10 kpc (|Haves et al.ll20l4 
lOstlin et al.ll20lll . Based on their curves of growth, a 
3-4 kpc diameter aperture would capture only one third 
to one half of the Lya flux. Indeed, aperture losses are 
confirmed for three galaxies in our sample. One Green 
Pea, 0926+4427 is also identified as LARS 14, and two 
(1133+6514 and 1219+1526) are serendipitously covered 
by GALEX grism observations. For the former, the COS 
aperture captured 40% of the total large aperture lumi¬ 
nosity estimated from LARS imaging ~ 1.4 x 10^^ 

erg s“^). The latter two cases show that 60 and 75% 
of the GALEX grism flux is detected by COS. Because 
the physical extent of the Lya emission may vary from 
galaxy to galaxy, and 7/10 of the Green Peas have no 
large-aperture measurements, we do not correct the COS 
measurements. Ultimately, though, the COS measure¬ 
ments are interesting because they tell us about the Lya 


that is observed within the central few kpc. These quan¬ 
tities are important for comparing to high-redshift spec¬ 
troscopic studies, where slits may subtend physical scales 
similar to the COS aperture. 

Besides the line fluxes. Table [3] lists the Lya equiv¬ 
alent width, luminosity, and the Lya escape fraction. 

The latter quantity is defined as the ratio of the 
observed Lya luminos ity to the intrinsic Lya luminosity, 
L°]^ya/L'^hya (s®® also iHaves et al.l [20141 1 . The intrinsic 
Lya luminosity is inferred from the dust-corrected Ha lu¬ 
minosity times the intrinsic case-B ratio of L^ya/Lua 
8.7[3 We also include measurements of the red and blue- 
side WLya for the 9/10 double peaked Lya lines. These 
quantities are calculated by directly integrating the emis¬ 
sion profiles for > 0 and < 0 separately. In two cases 
(0303-0759 and 1054+5238) the equivalent width of the 
blue peak is consistent with zero; this finding is not an 
indication that the blue peak is undetected. Rather, the 
blue emission is weak, and there is net absorption around 
zero velocity. 

Einally, Table[3]also lists kinematic signatures from the 
Lya lines. The velocities of the blue and red peaks are 
given, and the broad wings of the lines are quantified by 
calculating the velocity where Fx reaches the continuum. 

5. GOMPARISON TO NEARBY Lya SAMPLES 

The intrinsic Lj^y ^/ratio predicted by case-B recombi- 
na tion theory is often repor ted as 8.7 and incorrec tly attributed 
to IHumm^ Stor^ 1119871') or IBrocklehurstI 111971]) . Instead, the 
canonical assumption where 2/3 of all ionizing photons lead to a 
Lya photon yields Lya /H/3 = 23.1, and when Ha/H/3= 2.86 (for 
Te = 10, OOOi^), we have Lya/Ha = 8.1. However, the Lya /H/3 
ratio is density dependent, and 2/3 of recombinations lead to Lya 
photons only in the low density limit. At higher densities, collisions 
bump electrons from 2 ^5 to 2 ^P, thereby suppressing 2-photon 
contin uum and enhancing Lya emission. IDooita &: Sutherland! 
II20Q3I) tabulate Lya/H/5 and Ha/H/3 ratios for a range of tempera¬ 
tures and densities. For the present sample, the [S II] AA6716, 6731 
ratio indicates electron densities, rie, from 100 to 800 cm“^, im¬ 
plying intrinsic Lya/Ha=8.1 - 9.2. We adopt 8.7 as intermediate 
to these extremes. 


























In this section, we explore the role of dust and stars 
in regulating the Lya output of the Green Peas. In 
this context, it is useful to compare the Green Peas 
to nearby galaxies that have also been observed in 
Lya. We choose two samples with published Lya mea- 
surements: LARS (iHaves et al.ll2013L 12014 lOstlin et al.l 
I20l4 iRivera-Thorsen et all l2015h and nearby galaxies 
identified as Lya emitters (LAEs; W^ya > 15 — 20 
A) from GALE X Kr ism surveys (iDeharveng et al.l 120081 : 

iScarlata et al.l[20?M iGowie et al.ll2011h . We acknowledge 

that the selection effects and aperture size likely influence 
this comparison. On one hand, the GALEX LAEs were 
Lya selected via slitless grism spectroscopy with a large 
5" FWHM PSF. Under this configuration, we expect that 
most of the Lya emission should be included in the mea¬ 
surements, but the sample is biased towards rare, strong 
Lya emitting galaxies. On the other hand, the measure¬ 
ments reported for the Ha-selected LARS galaxies are in¬ 
tegrated inside Petrosian apertures that are defined using 
the Lya plus continuum images. These apertures, which 
range in diameter from 2.6 - 32 kpc, do not capture all 
of the Lya emission detected at large radii in the LARS 
images. For the low-mass, low-dust, high Wna end of 
the LARS sample, the typical 4 kpc diameter aperture 
is not too different than the COS aperture at z ^ 0.2. 
Indeed, for LARS 14/0926-1-4427, the Petrosian aperture 
Lya flux from the LARS image is comparable to (83% 
of) the the flux included in the COS spectrum. However, 
as we noted in (21 these measurements account for less 
than half of the total Lya flux in the LARS image. Nev¬ 
ertheless, we conclude that the LARS Petrosian aperture 
measurements and COS/Green Pea spectra are still use¬ 
ful, because they give a sense of how much Lya is escap¬ 
ing from the central regions of the galaxies. Finally, we 
acknowledge that Lya from ne arby COS-observed g alax- 
ies has also been presented by I Wofford et al.l (l2013fl , but 
we exclude this sample from comparison because their 
close proximity (z 0.03) implies that the COS aper¬ 
ture subtends only 0.6 kpc (unvignetted). Indeed, their 
COS NUV continuum images show much greater spatial 
extent compared to the Green Peas, suggesting that the 
Lya measurements are not easily compared to the other 
samples considered here. 

A straightforward approach for quantifying Lya escape 
is to directly compare Lya flux to the Ha and Ytp emis¬ 
sion that probe the nebular gas and dust. In Figure[31 we 
plot the Lya to Ha flux ratio against the Ha to H/3 flux 
ratio for the Green Peas, LARS galaxies, and GALEX 
LAEs. T he latter samp l e is c omp rised of 45 galaxie s 
taken from iScarlata et at] (1200911 and iCowie et al.l (1201111 . 
where slit-loss corrections facilitate the comparison be¬ 
tween Lya and Ha fluxes. Additionally, we show predic¬ 
tions from case- B recombination theo ry, as suming extinc¬ 
tion l aws from ICardelli et al.l (|1989ll and ICalzetti et al.l 
(|2000r) . In the absence of resonant scattering, the flux ra¬ 
tios should foll ow these lines. We als o show the clumpy 
dust models of IScarlata et al.l (|2009ll which form loops 
when the limit of high clump optical depth returns the 
observed line ratios to the intrinsic case-B ratios. While 
these models are needed to explain some of the high Lya 
/Ha ratios in dusty galaxies, they do not explain line 
ratios of the Green Peas which have low Lya to Ha for 
their dust content. Although the Green Peas show little 



Fig. 3.— The low Lya to Ha flux ratios of the Green Peas cannot 
be explained by dust. Solid and dashed lines show how, in the 
absence of resonant scattering, the intri nsic Lya to Ha flux ratio 
would respond to du st extinction for the ICalzetti et al.l <200(11 and 
ICardelli et al.l < 198911 extinc tion c urves. Additionally, clumpy dust 
models from IScarlata et al.l <2009l ~) are shown for 3, 5 and 7 clumps; 
they form loops when high-clump optical depth returns the flux 
ratios to their intrinsic, dust-free limits. Statistical errors for the 
Green Peas are small compared to the size of the data points, and 
are not reported for the other samples. Comparison with the LARS 
galaxies and the GALEX LAEs shows that the Green Peas have 
lower Ha/H/3 flux ratios, indicating less dust obscuration. GALEX 
LAEs identified as AGN have been excluded. 


to no dust, their Lya to Ha flux ratios span a factor of 
10. Two of the Green Peas, 1249-I-I234 and 09II-I-183I, 
fall close to the Lya to Ha ratio that is predicted for 
their dust extinction, but the remaining eight show ra¬ 
tios that are too low to be explained by dust extinction 
alone. 

To clarify how the Green Peas sample a different range 
of galaxy properties than LARS and the GALEX LAEs, 
Figure El shows histograms in stellar mass, SFR, Myy, 
UV-slope, (3, Wna, E{B - V)gas, LLya, the UV isopho- 
tal Petrosian radius, Ruv, as well as WLya, and 
The derivation of t hese properti e s are outlin ed in (J^lfor 
the Green Peas and IHaves et al.l (l20I3l . l2M3l for LARS. 
For the GALEX LAEs, we take the Muv, Who, LLya 
and WLya from Table I and 2 in Cowie et al. (2011; 
44 galaxies). The SFR and E{B — V)gas are calculated 
from the slit-loss co rrect ed line fluxes note d above (from 
IScarlata et al.ll2009l andiCowie et aI1l20Ill : 45 galaxies), 
assuming a iGalzetti et al.l (|200ClD exti nction law, an in - 
trinsic Ha/H/3 ratio of 2.86, and the iKennicutu (|1998r ) 
SFR calibration (converted to a iChabried 120031 IMF). 
GALEX LAEs identified as AGN have been excluded. 

Figure in shows that the Green Peas occupy a region of 
parameter space that is poorly sampled by other studies. 
Their stellar masses overlap with the low-mass end of the 
LARS galaxies, but extend to masses an order of magni¬ 
tude smaller. At the same time, the Green Peas’ SFRs 
and UV luminosities are, on average, higher than the 
LARS galaxies and GALEX LAEs. They are uniformly 
low in dust, with E{B — V)gas < 0.2 and P ^ —2.0. Fi¬ 
nally, Figure m shows that the sizes of the Green Peas are 
similar to the more compact half of the LARS galaxies. 
These different Green Pea properties seem to impact the 
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Fig. 4. — The Green Peas are compared to the LARS and GALEX 
LAE samples through their distributions of dust, stellar population 
properties, and Lya luminosity, equi valent width , and escape frac¬ 
tion. The LARS data are taken from IHaves et akl II20131 . 120141) . and 
the GALEX data are taken from Scarlata et al. (2009; SFRs and 
E[B — 1 ^) 503 ) and Gowie et al. (2011; Mjjv^ ^Hai and and 

Wli/q:)- The Ruv is the same isophotal Petrosian radius defined in 
IHaves et al.l II20131 . 120141) . measured in the FUV for LARS and the 
NUV acquisition images for the Green Peas. The distributions are 
normalized since the samples sizes being compared are varied. To 
our knowledge, stellar mass, /3 slopes, and UV sizes of the GALEX 
LAEs have not been published. 


Lya output: the Lya luminosities are an order of magni¬ 
tude larger than are observed for LARS and the GALEX 
LAEs. In fact, among these three samples, only the 
Green Peas have Lya luminosities in the r ange of most 
high redshift LAEs: L^ya > erg s“^ (|C)uchi et al.l 

[Moll . Moreover, WLya and show a broad range 
of values for the Green Peas, while the LARS galaxies 
are more peaked at low values. (The GALEX LAEs, by 
definition, exclude the low values otWLya and la 

fact. Figure [5] shows that in the combined set of LARS 
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Fig. 5. — The Lya escape fraction is strongly correlated with 
Wjyo, in the combined sample of Green Peas and LARS galaxies. 
The Spearman rank correlation coefficient is 0.64, and the prob¬ 
ability of the null hypothesis is 8.1 X 10”“*. The GALEX LAEs, 
on the other hand, show no significant correla tion. The LARS 
measurements are taken from IHaves et aLl II2014I L and the quanti- 
ties for the GALEX LAE s are calcula t ed fro m data tabulated in 
IScarlata et al.l Il2009l ~) and ICowie et al.l II2011I 1 . 


and the Green Peas, WHa is strongly correlated with 

Lya- The Spearman rank correlation coefficieniF^ is 
0.64 and a spurious correlation is rejected with a proba¬ 
bility of 8.1 X 10“^. 

This observation that the low masses, high SFRs, and 
low dust content of the Green Peas may fav or Lya emis¬ 
sion i s consistent with trends reported in I Gowie et al.1 
(j2011h . These authors compared GALEX LAEs to a UV- 
continuum selected control sample and found that the 
LAEs had bluer colors, more compact sizes, lower metal- 
licities, and higher Wua than their non-emitting coun¬ 
terparts. Similarly, Gowie et al. find Lya emission to 
be more common in samples with higher Wna-, and also 
reported a weak correlation between Wna and WLya- 
Here, by actually comparing /Ac“ instead of WLya for 
the optical emission line selected samples, we see that 
this correlation probably originates from increased Lya 
escape, rather than young stellar populations with intrin¬ 
sically high WLya- At the same time, however, we detect 
no statistically significant correlation between /Ac“ and 
WHa for the GALEX LAEs. 

Regardless, a relation between /Ac“ and Wna (or any 
of the other properties highlighted in FigureH]), does not 
explain the physical mechanism regulating Lya escape. 
Rather, it suggests that gas properties associated with 
the youth of a stellar population increase Lya escape. In 
the sections that follow, we use the ultraviolet interstellar 
absorption lines to further investigate the role of this gas. 

6. Lya EMISSION AND THE ABSORPTION STRENGTH OF 
INTERSTELLAR GAS 


The upper limits on (corresponding to net Lya ab¬ 

sorbers in the LARS galaxies) are set to zero in this statistical test. 
The Spearman correlation coefficient, using ranks, is not sensitive 
to upper limits in this case. 
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TABLE 4 

Observed UV Absorption and Emission Lines 


Ion 

F.' 

■‘-'ion 

(eV) 

Vacuum Wavelength 

(A) 

flu 

^ul 

(s-1) 

^low ^up 

(eV) 

HI Lye 

0 

937.80 

0.0078 

1.64 X 10® 

0.00 - 13.22070331 

H I Ly(5 


949.74 

0.014 

4.12 X 10® 

0.00 - 13.0545011 

H I Lyy 


972.54 

0.029 

1.28 X 10'^ 

0.00 - 12.7485388 

HI Ly/3 


1025.72 

0.08 

5.57 X 10'^ 

0.00 - 0 12.0875046 

H I Lyo 


1215.67 

0.46 

4.69 X 10® 

0.00 - 10.1988353 

O I 

0 

1302.17 

0.05 

3.41 X 10® 

0.00 - 9.5213634 

O I* 


1304.86 

0.05 

2.03 X 10® 

0.0196224 - 9.5213634 



1306.03 

0.05 

6.76 X 10® 

0.0281416 - 9.5213634 

Si II 

8.15 

1190.42 

0.277 

6.53 X 10® 

0.00 - 10.415200 



1193.29 

0.575 

2.69 X 10® 

0.00 -10.390117 



1260.42 

1.22 

2.57 X 10® 

0.00 - 9.836720 



1304.37 

0.09 

3.64 X 10® 

0.00 - 9.505292 

Si II* 


1194.50 

0.737 

3.45 X 10® 

0.035613 - 10.415200 



1197.39 

0.150 

1.40 X 10® 

0.035613 - 10.390117 



1264.74 

1.09 

3.04 X 10® 

0.035613 - 9.838768 



1309.28 

0.08 

6.23 X 10® 

0.035613 - 9.505292 

CII 

11.26 

1334.53 

0.129 

2.42 X 10® 

0.00 - 9.290460 

C II* 


1335.71 

0.115 

2.88 X 10® 

0.007863 - 9.290148 

Si III 

16.34 

1206.50 

1.67 

2.55 X 10® 

0.00 - 10.276357 

Si IV 

33.49 

1393.76 

0.513 

8.80 X 10® 

0.00 - 8.895697 



1402.77 

0.255 

8.63 X 10® 

0.00 - 8.838528 


Note. — Atomic line data are given for the transitions considered in this paper. 
Values are taken from the NIST Atomic Spectra Database. 


The most common approach to studying the impact of 
outflows on Lyo is to observe low ionization interstellar 
(LIS) metal lines. With ionization potentials less than 
13.6 eV, ions such as Si II and C II can trace both neu¬ 
tral and ionized hydrogen. As such, these transitions 
are typically used to quantify the H I gas that scat¬ 
ters Lya photons. In Figures [U [7] and [H we show the 
COS spectra covering Si II AA1190,1193, A1260, and C II 
A1334. The spectra have been normalized by linear fits 
to the local continuum, and zero velocity is marked by 
a dashed vertical line. Similarly, dotted lines mark the 
expected locations of fluorescent fine-structure emission 
lines, Si II* and C II*. These lines form when electrons 
excited by absorption in the Si II AA1190,1193, A1260, 
and C II A1334 transitions subsequently decay to the ex¬ 
cited ground state (followed by a fine structure transition 
to the ground state). We exclude from our analysis the 
O I and Si II lines at A1302, A1304, because their small 
wavelength separation and contribution from fluorescent 
and resonant re-emission (O I* A1304) complicates the 
interpretation (see spectra shown in the Appendix). We 
list the rest-frame vacuum wavelengths, ionization poten¬ 
tials, absorption oscillator strengths, /;„, and emission 
coefficients, A„/ for these lines in Tabled Additionally, 
we note that contamination to these ISM features by stel¬ 
lar absorption is negligible. In the top left panel of [6l [7] 
and [SI we show a 50 Myr old, c ontinuous star formin g 
model spectrum, with Z = 0.002 (jLeitherer et al.ir2014D . 
Similar stellar absorption is seen across a wide range of 
young, metal-poor stellar population properties, where 
UV spectra are 0-star dominated. 

Weaker equivalent width of these LIS absorption lines 
have been associated with increas ed Wr,„rv in stacked 
spectra of galaxi es at z ^ 3 — 4 (iShaolev et al.l 120031 : 
[Jones et al.l 1201211 . Qualitatively, the Green Peas lend 
some support to this scenario: the two strongest ab¬ 
sorbers 1244-1-0216 and 1137-1-3524 are among the lower 


/e^c“ Green Peas, and the weakest absorbers (1133-1-6514 
and 1219-1-1526) are among those with the highest 
To quantify this trend, we measure equivalent widths of 
the LIS metal absorption lines by directly integrating the 
normalized spectra out to the velocity where the absorp¬ 
tion meets the continuum. The equivalent width uncer¬ 
tainties are determined by propagating the error vector 
over the same velocity range, including a systematic 10% 
uncertainty on the continuum normalization. For unde¬ 
tected lines, we approximate the upper limit by taking 
the equivalent width of the marginally detected lines. For 
the LIS lines this threshold is around 0.5 (0.4) A in the 
observed (rest) frame. These measured equivalent widths 
are listed in Tabled 

In FigurelH we compare W hya to the equivalent widths 
of Si II A1260 and G II A1334 for both t he Green Peas 
and th e stacked z ~ 3 LEG sample from IShaolev et ahl 
(|2003f) . In this plot, the Green Peas are consistent with 
the z ~ 3 LBGs, although they extend to higher W hya 
and weaker Wlis- Nevertheless, we do not detect a sig¬ 
nificant correlation; the Spearman correlation coefficient 
(p) and probability of the null hypothesis (P) are given in 
the upper left of each panel in Figure [5] However, since 
'WLya is only a rough proxy for Lya escape, we recast 
this diagram using in the right panels of Figure O 
Here, the significance of the correlation is improved: the 
Spearman coefficient (which, using ranks, is not sensitive 
to the non-detections in this case) implies that the corre¬ 
lation is robust at 91% confidence for Si II and 97% confi¬ 
dence for G II. Although these trends are only marginally 
significant, the consistency between G II and Si II sup¬ 
ports a real correlation. Moreover, these trends work in 
the same direction as the z 3 — 4 measurements from 
stacked spectra. It appears that conditions which create 
weaker LIS absorption lines also favor greater /^^“. We 
will return to this trend and discuss its possible meaning 
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Fig. 6 .— The COS spectra of 10 Green Peas show weak, 
blueshifted absorption in Si II 1190.4 and 1193.3. The velocity 
scale is appropriate for Si II A 1193.3, with zero velocity marked 
by the vertical dashed line. The vertical dashed line around -700 
km s“^ shows the expected velocity for absorption in Si II 1190.4, 
and the Si II* emission at 1194.5 A and 1197.4 A is marked by 
dotted lines. The portions of the spectra that are plotted in grey 
are affected by Milky Way absorption or extra noise, except for the 
case of 1133+6514 where we suspect intervening Si II A 1260 from 
an absorber at 2 = 0.17. The horizontal dashed line shows unity in 
the normalized spectra. The blue spectrum, overplotted in the top 
left panel, is a Z = 0.002 Starburst99 model, which confirms that 
stellar Si II absorption is absent from young, UV-luminous stellar 
populations. The spectra are displayed in order of increasing 
(left to right, top to bottom). 

in where we offer a more comprehensive interpreta¬ 
tion of the data. 

Despite heavy reliance on the LIS metal lines, they 
remain indirect probes of the neutral hydrogen that is 
scattering the Lya emission. Fortunately, the COS spec¬ 
tra of the Green Peas allow a more direct look at the H I 
gas, with at least one Lyman series line (besides Lya) 
observed in each of the 10 galaxies. Figure ITUl highlights 
these features, showing either Ly/3 or Lyy for each Green 
Pea. Whenever more than one Lyman series line is ob¬ 
served (7/10 galaxies) we find that their absorption pro¬ 
files are consistent. The equivalent widths of these lines 
(Table O are measured in the same manner as for the 
metal LIS lines, although we instead adopt a 20% un¬ 
certainty on the continuum placement because normal¬ 
ization is more challenging at these wavelengths. Addi¬ 



Fig. 7.— Same as Figure [B] but for Si II A 1260.4 and Si II* 
A 1265.0. These lines fall beyond our wavelength coverage for 
1424+4217. 

tionally, the same stellar absorption model used for the 
metal lines is shown for Ly/3 and Lyy. In contrast to the 
metals, the stellar H I absorption is significant around 
zero velocity. Nevertheless, the stellar absorption does 
not explain the blueshifted and highly opaque absorp¬ 
tion, indicating a significant contribution from outflow¬ 
ing interstellar gas. We conclude that these lines are still 
useful for probing the outflowing gas at moderate to high 
(blueshifted) velocities. 

The H I absorption lines show some remarkable dif¬ 
ferences from the LIS metal lines. First, the equivalent 
widths of the Lyman series lines show no variation with 
^Lya- Second, and in contrast to the spec¬ 
tra shown in Figures HI [7] and HI the hydrogen absorp¬ 
tion lines show little to no residual intensity at mod¬ 
est blueshifted velocities (a few hundred km s“^). This 
finding indicates that near 100% of the stellar light is 
covered by H I absorbing gas at these velocities. More¬ 
over, the high opacity and the similarity between Ly/3, 
Lyq, hyS and Lye equivalent widths (when more than 
one line is detected) indicates that these lines are satu¬ 
rated but not damped. Hence, the neutral hydrogen col¬ 
umn density is poorly constrained: Nhi ^ 10^® — I0^° 
cm“^. Remarkably, even the low end of this range leads 
to high Lya optical depth at line center, even if it is op¬ 
tically thin to hydrogen-ionizing Lyman continuum pho- 
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Fig. 8. — Same as Figures [0 and [71 but for C II A 1334.5 and C II* 
A 1335.7. Again, these lines fall beyond our wavelength coverage 
for 1424+4217. For 1054+5238, the C II A 1334.5 absorption is 
heavily contaminated by strong Milky Way A1 II A1671. 


tons > 10^ — 10^: iVerhamme et al.ll2006L I2014D . 

Finally, we show higher ionization Si III A1206.5 and 
Si IV A1394, A1403 in Figures ITT] and IT^ The equivalent 
widths of these lines are also listed in Table [5] Unlike the 
LIS metal lines, the equivalent widths of these higher ion¬ 
ization lines do not appear to change across the sample. 
This conclusion is most apparent for the Si III lines, all 
of which have good S/N and IT ~ —1.3 A. In this as¬ 
pect, t he Green Peas are al so consistent with the z ~ 3 
LBGs: iShaplev et al.l (l2003f) report that, despite the re¬ 
lation between WLya and Wlis^ Wsi iv is the same in 
each of their four WLya-defined stacks. 

7. DO OUTFLOWS HELP THE Lya PHOTONS ESCAPE? 

The kinematics of the CGM may play an important 
role in regulating Lya escape from galaxies. When the 
photons scatter in outflowing gas, they can Doppler shift 
out of resonance with the bulk of the ISM and escape 
more easily. Hence, we may expect that a large veloc¬ 
ity gradient in the CGM promotes strong Lya emission. 
Early studies find evidence for this scenario. In a sam - 
ple of eight local starburst galaxies, iKunth et al.l (|1998[1 
find that half of their sample with damped Lya absorp¬ 
tion showed no evidence of outflowing gas in their LIS 
lines. Yet the other half of their sample that showed 


Lya emission exhibite d outflows with veloc ities around 
200 km s“^. Similarly. iWofford et al.l (|2013ll report sim¬ 
ilar findings from COS observations of 20 nearby galax¬ 
ies. In their sample, seven Lya emitting galaxies have 
vlis ~ —100 km s“^, whereas ten galaxies with damped 
Lya p rofiles have t’r.r.s' —20 km s“^. And most re¬ 
cently, iMartin et al.l ( 201511 show that Lya escape from 
ULIRGs is enhanced when the wings of the Lya and 
[O III] lines reach greater blueshifted velocities. 

In this section, we will test whether the Lya escape is 
aided by outflowing gas in the Green Peas. To quantify 
the outflow kinematics, we will use metal and hydrogen 
absorption lines, and also explore the kinematic measures 
from the Lya emission. 


7.1. Gas Kinematics in Absorption 

Figures m Hand m illustrate the kinematics of the LIS 
ions. When the lines are detected, they are blueshifted 
with wings extending several hundred km s“^. To mea¬ 
sure the gas kinematics, we estimate characteristic out¬ 
flow velocities, Vc- Since many of the lines have non- 
Gaussian shapes, we avoid fitting centroid velocities, and 
instead calculate the equivalent width weighted velocity: 


Vc = 


J t(1 Jnorrn)d\ 
/(I - fnorm)dX 


( 1 ) 


where fnorm is the normalized flux, v represents velocity 
at each wavelength, and the denominator is easily recog¬ 
nized as the equivalent width. The resultant character¬ 
istic velocities are listed in Table [T] The errors on these 
velocities are calculated using 1000 Monte Carlo realiza¬ 
tions where the spectrum is perturbed according to its 
error vector and a 10% continuum normalization uncer¬ 
tainty. We also quantify the maximum outflow velocity 
by determining where the absorption trough reaches the 
continuum. These velocities are listed in Table [SJ the 
uncertainties are calculated using a Monte Garlo simula¬ 
tion in the same manner as for the characteristic outflow 
velocities. For completeness, and to obtain some kine¬ 
matic information for 1133-1-6514 and 1219-1-1526 which 
have no LIS lines, we also list the characteristic and max¬ 
imum outflow velocities for the higher ionization Si III 
and Si IV lines. And likewise, maximum outflow veloci¬ 
ties for the Lyman series lines are listed in Table |9] (where 
we have again included a 20% continuum normalization 
uncertainty). However, because of the stellar component 
and the breadth of the H I absorption lines, we do not 
measure characteristic outflow velocities for these fea¬ 
tures. 

Figure [13] shows how depends on these kinematic 
features for the H I, LIS metals, and higher ionization 
lines. The top two panels focus on Vc, with different 
symbols indicating different transitions. The upper left 
panel shows that we detect no significant trend between 
and outflow velocity measured from the LIS met¬ 
als. However, we must note that two of the three high¬ 
est galaxies cannot be included in this plot and 

correlation test, because they are undetected in all LIS 
absorption lines (1133-1-6514 and 1219-1-1526). Their ad¬ 
dition in the upper right panel of Figure 1131 where we 
show the higher ionization states of Si, hints at a cor¬ 
relation. For Si III, which is detected at good S/N in 
all ten Green Peas, the Spearman correlation coefficient 
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^Si II A1260 




II A1334 II A1334 (A) 


Fig. 9.— The rest-frame equivalent widths of Lyo; emission (left panels) and the Lya escape fraction, feJi^ (right panels) are plotted 
against the absorption equivalent width of the low ionization Si II A 1260 and C II A1334 lines. The Spearman correlation coefficient, p, 
and probability of the null hypothesis, P, are giv en in the upper left of each panel. Comparison data from stacked spectra of 2 ~ 3 galaxies 
are shown as grey squares IlSnaplev et al.ir20Q^3l) . While the Green Pea data are consistent with an extrapolation of the 2 ~ 3 data, this 
correlation is only significant when feJc^^ is used to quantify Lya output. 


rules out a spurious relation at 98% confidence. But this 
trend is only significant with the addition of the two high 
/isc“ Green Peas. In the bottom panels, where we show 
the maximal velocities, the metals tell a similar story: 
the Green Pea 1219+1526 (the highest object) has 
greater maximal outflow velocity than most when Si III is 
considered. But otherwise, we detect no trend. Finally, 
in the bottom left panel of Figure [U] we add maximal 
outflow velocities reached by H I Lyman series lines (red 
points). These lines consistently reach —700 to —800 km 
s“^, and show no trend with In summary, since 

the only possible correlation detected in Figure [T3] re¬ 
lates with higher ionization lines that do not trace 
H I, we conclude that these data give no compelling ev¬ 
idence for a scenario where Lya escape is enhanced by 
scattering in outflowing H I gas. 


Remarkably, Figure [13] shows that the maximal out¬ 
flow velocities measured by the H I Lyman series lines 
are systematically bluer than the maximal velocities of 
the LIS metal lines. This result demonstrates that the 
H I absorption is markedly more sensitive to low den¬ 
sity, high velocity gas that cannot be detected in the LIS 
metal absorption lines. Yet this H I likely plays an im¬ 
portant role in scattering Lya photons; we will explore 
this topic more in ED where we compare the Lyman 
series absorption and Lya emission kinematics. 

7.2. The velocity structure of Lya emission 

The Lya emission line profiles give another probe of the 
outflowing gas kinematics, as the broad emission is (at 
least in part) generated by resonant scattering in the gas 
around the galaxy. Table [3] gives kinematic measures of 
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TABLE 5 

Equivalent widths of ISM metal absorption lines 


Object ID 

Si II A1190 

Si II A1193 

Si II A1260 

C 11 A1334 
(A) 

Si III A1206 

Si IV A1393 

Si IV A1403 

0303-0759 

-0.4 ±0.3 

-0.4 ±0.3 

-0.6 ±0.2 

-1.2 ±0.3 

-1.3 ±0.2 

-1.1 ±0.4 

-0.9 ±0.4 

1244+0216 

-1.1 ± 0.3 

-1.0 ±0.3 

-1.6 ±0.4 

-1.3 ±0.5 

-1.1 ±0.3 

-1.3 ±0.8 

-1.0 ±0.8 

1054+5238 

-0.8 ±0.3 

-0.5 ±0.3 

-1.0 ±0.5 


-1.8 ±0.5 

-1.9 ±1.2 

-1.3 ± 1.2 

1137+3524 

-1.0 ±0.3 

-1.1±0.3 

-1.2 ±0.3 

-1.3 ±0.3 

-1.4 ±0.2 

-1.3 ±0.3 

-1.5 ±0.5 

0911+1831 

-0.6 ± 0.2 

-0.9 ±0.2 

-1.1 ± 0.4 

-0.7 ±0.4 

-1.3 ±0.4 

-1.9 ±0.7 

> -0.6 

0926+4427 

-0.4 ±0.1 

-0.3 ±0.1 

-0.5 ±0.1 

-0.7 ±0.2 

-1.5 ±0.2 

-1.7±0.3 

-1.3 ±0.2 

1424+4217 

-0.7 ±0.4 

-0.3 ±0.3 



-1.0 ±0.3 



1133+6514 


> -0.4 

> -0.4 

> -0.4 

-1.3 ±0.4 

> -0.6 

-0.6 

1249+1234 

-1.8 ±0.4=- 

-0.9 ±0.3 

-0.8 ±0.4 

-1.1 ±0.2 

-1.3 ±0.5 

-0.9 ±0.6 

1219+1526 

> -0.4 

> -0.4 

> -0.4 

> -0.4 

-1.2 ±0.3 

-1.6 ±0.6 

-0.7 ±0.4 


Note. — Rest-frame equivalent widths of metal absorption lines discussed in this paper. 
^ This equivalent width is the total for the blended Si II A1190, A1193 lines. 


TABLE 6 

Equivalent widths of H I absorption lines 


Object ID 

Ly^ 

Ly7 

(A) 

Ly5 

Lye 

0303-0759 

-2.6 ±0.6 




1244+0216 

-2.6 ± 0.7 


-2.1 ±0.6 


1054+5238 

-2.8 ± 0.7 


-2.5 ±0.6 

1.9 ±0.5 

1137+3524 

-2.5 ± 0.8 

-2.5 ±0.6 



0911+1831 

-1.8 ±0.7 



1.2 ±0.5 

0926+4427 

-2.4 ±0.5 




1424+4217 

1133+6514 

-2.9 ± 0.7 

-1.9 ±0.5 

-2.0 ±0.7 


1249+1234 

-2.2 ± 0.7 

-2.1 ±0.6 


2.2 ±0.8 

1219+1526 

-2.0 ±0.7 

-2.6 ±0.7 




Note. — Rest-frame equivalent widths of H I Lyman series 
absorption lines discussed in this paper. Lyman series equivalent 
widths contain a contribution of around 0.6 to 0.8 A from stellar 
absorption. 


TABLE 7 

Characteristic outflow velocities from ISM absorption lines 


Object ID 

?^1190 

L'1193 

?^1260 

V1334 

L’1206 

1^1393 

V1403 

0303-0759 

a 

a 

-180 + 60 

-220 ± 60 

-240 ± 50 

-240 ± 60 

-230 ± 80 

1244+0216 

-100 + 20 

-50 ± 20 

-110 + 20 

-80 ± 20 

-70 ± 20 

-230 ± 70 

-100 + 50 

1054+5238 

-140 + 100 

-120 + 80 

-280 ± 70 

b 

-200 ± 40 

-180 + 80 

-160 ± 70 

1137+3524 

-140 + 30 

-130 + 30 

-190 + 30 

-150 + 40 

-170 + 40 

-140 + 30 

-200 ± 40 

0911+1831 

-200 ± 80 

-300 ± 90 

-360 ± 70 

-130 + 90 

-400 ± 90 

-290 ± 50 

-140 ± 100 

0926+4427 

a 

a 

-330 ± 170 

-230 ± 100 

-320 ± 60 

-320 ± 70 

-270 ± 70 

1424+4217 

-260 ± 80 

-210 + 90 

C 

C 

-280 ± 60 

C 

C 

1133+6514 

a 

a 

a 

a 

-270 ± 50® 

a 

a 

1249+1234 

d 

d 

-150 + 60 

-210 + 60 

-250 ± 60 

-210 + 50 

-250 ± 80 

1219+1526 

a 

a 

a 

a 

-420 ± 90 

-370 ± 80 

a 


Note. — Equivalent width weighted velocities are given for Si II AA1190.4,1193.3, A 1260.4, C II A 1334.5, 
Si III A 1206.5, and Si IV AA 1393.8, 1402.7. 


^ Line is marginally detected or completely undetected, so reliable kinematic information is unavailable. 

^ Line is contaminated by Milky Way absorption. 

^ Line is not covered by the present spectra. 

Si II AA 1190.4, 1193.3 lines are blended. 

® The blue wing of the line may exten d m ore than —1200 km s~^, but uncertainties in continuum normalization 
make this extent unclear (see Figure EUi. We report the velocity from a Gaussian fit to the main absorption 
component of this line. 
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TABLE 8 

Maximum Velocities of Metal Absorption Lines 


Object ID 

„.max 

M190 

..max 

M193 

'1^1260 

..max 

^1334 

..max 

^1206 

..max 

^1393 

..max 

^1403 

0303-0759 

a 

a 

b 

-560 + 50 

-710 + 100 

-580 + 120 

-440 + 30 

1244+0216 

-410 + 60 

-330 + 70 

-550 + 90 

-350 + 60 

-320 + 80 

-630 + 90 

-400 + 70 

1054+5238 

-440 ± 60 

-420 + 100 

-760 + 150 

b 

-690 + 100 

-650 + 110 

-490 + 70 

1137+3524 

-440 ± 90 

-400 + 60 

-590 + 130 

-450 + 40 

-510 + 110 

-500 + 100 

-630 + 160 

0911+1831 

-410 + 110 

-600 + 50 

-800 + 200 

-600 < i; < -300^ 

-930 + 240 

-700 + 160 


0926+4427 

a 

a 

-580 + 130 

-580 + 90 

-880 + 100 

-930 + 90 

-680 + 60 

1424+4217 

-630 + 130 

-300 + 70 

C 

C 

-560 + 70 

C 

C 

1133+6514 

a 

a 

a 

a 

-600=^ 

a 

a 

1249+1234 

-500 + 50 

e 

-390 + 70 

-420 + 90 

-540 + 50 

-470 + 120 

-400 + 110 

1219+1526 

a 

a 

a 

a 

-960 + 160 

-670 + 220 

a 


Note. — Maximum velocities (in km s“^) are given for the same lines as in Table[7| 

^ Line is marginally detected or completely undetected, so reliable kinematic information is unavailable. 

^ Line is contaminated by Milky Way absorption. 

^ Line is not covered by the present spectra. 

^ The blue wing of the lin e ma y extend more than —1200 km s“^, but uncertainties in continuum normalization make this 
extent unclear (see Figure FTT|l . Most of the absorption, however, lies within —600 km s“^. 

® The continuum is not reached before the absorption from neighboring Si II A1190.4 is reached. Therefore the maximum 
outflow velocity from this line is unreliable. 

^ This absorption profile appears divided into two parts, possibly from a noise spike that crosses unity in the normalized 
spectrum around -300 km s“^. The maximum velocity in the more blueshifted portion of the profile is around -600 km 
s“^. 


TABLE 9 

Maximum Velocities of Lyman Series Absorption Lines 


Object ID 

..max 

^LyP 

..max 

Lyj 

..max 

^LyS 

..max 

^Lye 

0303-0759 

-800 <v< -200=^ 




1244+0216 

-660 + 70 


-530 + 40 


1054+5238 

1137+3524 

-800 + 60 

-750 + 90 

-670 + 50 


0911+1831 

-760 + 180 

-1200 + 270 


-730 + 130 

0926+4427 

1424+4217 

-700 + lOO*" 

-700 + 110 



1133+6514 

-720 + 60 

-800 + 170 

-820 + 160 


1249+1234 

-620 + 100 

-690 + 90 


-760 + 170 

1219+1526 

-660 + 70 

-740 + 30 




Note. — Maximum velocities, (in km s“^) are given for the Lyman series 
lines. Lines which are not listed are either contaminated by Milky Way ab¬ 
sorption, fall beyond the blue-wavelength cutoff of our COS spectrum, or near 
the gap between the FUV segments, tlnlike the metal lines, we do not fail to 
detect absorption from H I gas. 

^ Blue wing of line impacted by Milly Way absorption. 

^ The normalized continuum does not reach unity before nearby Milky Way 
absorption sets in (see Figure fTOll . Visual inspection of the spectrum suggests 
an approximate maximum velocity of around —700 di 100 km s“^. 










16 



Fig. 10.— Lyman series absorption line profiles directly trace 
hydrogen gas. We show the Ly^ absorption line for 8/10 Green 
Peas, and Lyy for 2/10. The latter have Ly/3 that is contaminated 
by geocoronal emission (1424+4217), or fell at a noisy wavelength 
near the gap between the COS FUV segments (1249+1234). Star- 
burst99 model spectra, illustrating the stellar component of the 
Ly/3 and Lyy lines are shown in blue. The model plotted is a 50 
Myr, continuous star-forming population with Z = 0.002; simi¬ 
lar profiles are present across a wide range of stellar population 
models, provided that UV spectra are O-star dominated. As with 
Figures 0 [3 and [H] the grey shaded regions of the spectra show 
contamination from Milky Way absorption, geocoronal emission, 
and in the case of the Lyy profiles, the adjacent C III A977.0 line. 

the Lja lines, including the velocities of the red and blue 
peaks, and the maximal velocities of the red and blue 
wings of the lines. Figure [H] shows how changes 
with these Lya kinematic measures. Statistically signif¬ 
icant correlations appear in some, but not all, of these 
quantities. Each panel in Figure [TJ] lists the Spearman 
correlation coefficient, p, and the probability, P, that 
the correlation arises by chance. For this analysis, we 
exclude the Green Pea with single-peaked Lya emission 
(1249-1-1234) from the statistics and plots involving the 
blue peak and the peak separation, but it is included 
in the other diagnostics. Similarly, the maximal red ve¬ 
locity for 1424-1-4217 is a lower limit (the spectrum was 
truncated due to a failed observation), and this galaxy is 
excluded from the relevant statistics. 

Figure [HI shows that the separation between the Lya 
emission peaks becomes smaller when is larger (top 
right panel), and that this trend is driven primarily by 



a shift in the blue emission peak (top left panel). At 
the same time, the galaxies with higher may show 
somewhat lower velocities for their red Lya peaks, but 
the correlation is not significant. Compared to the blue 
emission peaks, the red Lya emission peaks also inhabit a 
smaller range of velocities. These findings are consistent 
with studies of high-redshift galaxies (z ^ 2 — 3), where 
the kinematics of the red Lya emissio n peak has been 
studie d by multiple group s. Not ably, iHashimoto et al.l 
lj2013ll and iShibuva et al.l (j2014ll measure a red peak 
velocity, AvLya ~ 200 km s“^ for 10-20 LAEs with 
^Lya> 50 A, compared with 450 km s“^ measured for 
UV-continuum sel ected Lyman Break Galaxies (LBGs; 
iSteidel et al.ll20l5l . In this sense, the Green Peas are 
more similar to the high-redshift LAEs than with the 
LBGs. Indeed , for a larger sample of 158 galaxies, 
lErb et akl (120141 1 report a significant anti-correlation be¬ 
tween Wiya and AvLya, such that Lya is stronger when 
it emerges closer to the systemic velocity. Our trend of 
increasing with peak separation fits naturally with 
this scenario; increased emission near the systemic veloc¬ 
ity serves to shift the emission peaks closer together. A 
visual inspection of Figure [2] supports this idea; where 
the Lya emission is strong, there is net emission around 
ti ~ 0, whereas net absorption appears around the sys¬ 
temic velocity when is low. Here, for the first time, 
we show that the kinematics of the blue Lya emission 


































































































17 



Fig. 12.— Same as[6l[7][8l and llll but for Si IV AA1393.8,1402.8. 
The velocity scale is appropriate for Si IV A1393.8, and the dashed 
line around 1900 km s“^ marks zero velocity for Si IV A1402.8. 


peak show even greater variation than the red peak. 

This tight correlation between and peak separa¬ 
tion suggests that Lya escape is determined by the neu¬ 
tral hydrogen column density. In short, when the column 
density is low, the Lya photons can escape nearer to the 
systemic velocity, with less scattering in the expanding 
circumgalactic envelope. Variations in dust, on the other 
hand, preserve the shape of the profiles while decreas¬ 
ing the strength of the emission (|Verhamme et al1l2006L 
120081 [201^ IBehrens et al.ll2014n . Indeed, the H I density 
has already been used to explain the trends with the red 


cussed above; 

Gtionis et al.ll2013l: iHashimoto et al.ll2013l: 

IShibuva et al. 

12014 |Erb et al.ll2014D. Moreover, recent 


show that small separations between peaks. An < 300 
km s“^, imply Nhi ^ 10^® cm“^, whereas higher sepa¬ 
rations 300 < An < 600 km s“^, arise in galaxies with 
10^® < < 10^0 cm“^. 

However, the Lya ra diative transfer models models of 
iVerhamme et al.l (120141) do not produce all the features in 
our data. First, unlike the Green Peas, increasing N^i in 
the models shifts the red peak to higher velocities while 
holding the blue peak at a nearly fixed velocity. Second, 
in modeled profiles, the high-velocity wings of the Lya 
reach only to ± a few hundred km s“^ in contrast to 


the several hundred km s“^ observed in the Green Peas. 
This discrepancy may be attributable to an assumed in¬ 
trinsic Lya profile which neglects the bro ad wings ob¬ 
served in the optical emission lines. Indeed. [Martin et al.l 
(j20I5D successfully produce broad Lya wings on eight 
ULIRGs by modeling the Lya emission as a superposi¬ 
tion of intrinsic (broad + narrow) and scattered compo¬ 
nents. And the third difference between the data and 
models of Verhamme et ah: the double-peaked emission 
profiles appear mostly when the shell expansion veloc¬ 
ity is low, V < 100 km s“^. Yet the Green Peas show 
—300 ^ Vc —100 in their LIS lines. Nevertheless, 
Lya radiative transfer mo deling offers some intrig uing 
lines of investigation. As IVerhamme et al.l (|2014 ) al¬ 
ready noted, the galaxies with small peak separations 
may have low enough column densities to be optically 
thin to the hydrogen-ionizing Lyman continuum (LyC) 
photons. The Green Peas with the most closely spaced 
peaks, 1219-1-1526 and 1133+6514, make excellent candi¬ 
dates for followup observations aimed at direct detection 
of the LyC. Remarkably, among the present sample of 
Green Peas, 1133+6514 shows lower Who and W Lya de¬ 


spite its high These conditions are also consistent 

with LyC leakage. 

In addition to the peaks of the Lya emission, we ex¬ 
plore the kinematics probed in the wings of the lines. 
The bottom panels of Figure [14] show how the maximal 
velocity reached in the Lya lines relates to In 

these diagnostics, the typical 200 km s“^ uncertainty on 
the maximal velocities makes it difficult to draw robust 
conclusions. Although only the center panel (illustrat¬ 
ing the maximal velocity in the red wing of the line) 
shows a significant relation, these plots hint at a sce¬ 
nario where broader Lya emission wings may be associ¬ 
ated with increased Lya escape. This finding is qual¬ 
itatively c onsistent with the sample of eight ULIRGs 
studied by iMartin et al.l (|2015D , where increased Lya es¬ 
cape was associated with greater blueshifted Lya and 
[O III] velocities. In this case, the close association of 
the blueshifted Lya and [O III] A5007 kinematics sug¬ 
gested an in-situ production of Lya photons in the high 
vel ocity gas (rather th an scattering). As demonstrated 
by IMartin et al.l (j2015f ). this emission would result from 
cooling of the hot galactic wind. However, with the reso¬ 
lution of the SDSS spectra, it remains difficult to test 
whether the Lya kinematics in the wings correspond 
closely with the [O III] and Ha kinematics. Ultimately, 
more work is needed to obtain higher spectral resolution 
observations of the nebular gas, so that we may compare 
the line profiles in greater detail. We defer this analysis 
to a future study. 

In order to more fully understand what the Lya kine¬ 
matics are measuring, we also compare the spectra to 
H I absorption profiles measured by the Ly/3 or Lyy lines 
in Figure 1101 This comparison shows that the maxi¬ 
mal velocities of the Lya and H I absorption are well 
matched (except for possibly 1244+0216), extending to 
around —700 km s“^ in most cases. The presence of 
H I absorption over the same range of velocities as the 
blueshifted Lya emission is noteworthy; it demonstrates 
clearly that the gas which scatters Lya photons exists in 
an envelope spanning a range of velocities. Shell models, 
where the cool gas exists at a single velocity, cannot de- 
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Fig. 13.— The Lyo: escape fraction shows no convincing evidence for a correlation with outflow velocity in the low-ionization metals 
or H I (left panels), but hints at a possible trend with the velocity of more highly ionized gas (right panels). The top panels show the 
equivalent width weighted characteristic outflow velocities of different ions (Table 0, with a typical 60 km s“^ error bar plotted in grey. 
The bottom panels show the maximal outflow velocity (Tables 0 and O, measured from the velocity where the absorption profile meets 
the continuum. The typical errors in this case are around 100 km s“^. The bottom left panel includes the H I absorption (red points) 
which occupy a consistently narrow range of maximal velocities (~ —600 to —800 km s“^), and are systematically bluer than the metal 
line tracers that are often use to probe H I. 


scribe absorption spanning several hundred km s“^, even 
though they often provide satisfacto ry fits to Lya emis- 
sion profiles (iVerhamme et alJl2008h Pre vious studies 
(e. g. lPettini et al.ll2002l : iKulas et al.ll2012ll have already 
noted this shortcoming of the shell model, drawing a sim¬ 
ilar conclusions from the velocity gradient seen in metal 
lines. Here, for the first time, we are able to demon¬ 
strate the need for an envelope of H I gas, without the 
need for metals as a proxy. This added constraint, it 
turns out, is important since the LIS and Lyman series 
lines show somewhat different kinematics in their maxi¬ 
mal outflow velocities. As we showed in Figure [131 the 
LIS metal lines reach maximal velocities around —400 to 
—600 km s“^, whereas the Lyman series lines extend to 
between —600 and —800 km s“^. We interpret this differ¬ 
ence as an indication that the LIS metal absorption lines 
are insensitive to the gas that constitutes the wings of 
the Lyman series lines. Nevertheless, the excellent corre¬ 


spondence between the blueshifted Lya emission and H I 
absorption velocities suggests that this high velocity, low 
density gas is still important for scattering Lya photons 
and creating the profiles that we observe. 

Stacking analyses of UV-luminous galaxies at z ^ 2 —4 
have also explored the ori gin of the blue peak emission by 
comp aring it to LIS lines ([Steidel et al.ll2010l : IJones et"^ 
1201211 . These studies have found that the velocity of the 
blue peak most closely corresponds with the velocities in 
the blue wing of the LIS absorption lines, around —500 
to —600 km s“^. In other words, the blue Lya emis¬ 
sion is strongest not where the apparent optical depth of 
the LIS lines is greatest, but instead where it is decreas¬ 
ing. The Green Peas show that these characteristics of 
stacked spectra is not uniform among galaxies. In Figure 
1161 we summarize the kinematic measurements from H I 
absorption, LIS metal absorption (Si II and C II), and 
blueshifted Lya emission. We show the velocity range 
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Fig. 14. — The velocity of the blue Lya peak, as well as the velocity separation between the blue and red peaks are strongly correlated 
with the escape of Lyo; photons. Other measures, including the maximal velocities, and the velocity of the red peak are only weakly related 
to Lya escape. Each panel is labelled with the Spearman correlation coefficient, p, and P, the probability that the correlation arose by 
chance. The grey error bars show the typical uncertainty in each measurement. 


probed by each set of features, by computing the error 
weighted mean and its uncertainty. The Green Peas are 
sorted by , with a comparison from st acked spectra 
of LB Gs at z 2—4 shown a t bottom ffrom ISteidel et ^ 
l201fll and I.Tones et'H] 12012(1 . While the LIS absorption 
velocities are similar between the high redshift samples 
and Green Peas, the latter show blueshifted Lya emis¬ 
sion emerging closer to the systemic velocity and Vc- In 
fact, while the lower Green Peas show the blue peak 
emission falling between the LIS Vc and Vmax, this trend 
may break down as increases. Both 0926-F4427 and 
I424-F4217 show that the blueshifted Lya peak probably 
emerges at lower blueshifted velocities than Vc- More¬ 
over, if the LIS metals which are too weak to detect in 
11334-6514 and 12194-1526 still follow the kinematics of 
the more highly ionized gas, we might expect this trend 
to be stronger. These galaxies show Vc (Si III) of —270 
and —420 km s“^ respectively, compared to their Lya 
■^peal; around -100 km s“^. 

In summary, we conclude that while the wings of the 
Lya emission lines emerge over the same range of veloc¬ 
ities probed by H I absorption, the peak emission shows 
little relation to any of the outflow kinematic measures. 
This finding is consistent with the conclusions from pre¬ 
vious high-redshift studies, where it is argued that the 
emergent Lya emission is more strongly affected by the 
optical depth of gas near the systemic velocity than by 
outflow kinematics (iSteidel et a.l.l[201Ct iLaw et aT]l2012l : 
iGhonis et al.ll201^ lErb et al.ll2014lb 


8 . SUMMARY & DISCUSSION 

In our analysis so far, we have explored the role of ISM 
and CGM gas in regulating the amount of Lya emission 
that we observe from the Green Peas. In this section 
we will discuss the physical implications of our measure¬ 
ments. First, however, we summarize the findings from 
the previous sections: 

1. The “Green Pea” classification identifies objects 
with prominent Lya emission and little to no Lya 
absorption. Nine of the ten galaxies in our sam¬ 
ple have Lya profiles that show both blue and red 
peaks, and unlike other low redshift samples, their 
Lya luminosity and W Lya reach the r ange probed 
by high redshift Lya suryey s (e.g. lOuchi et al.l 
lamil: [Henry et a1J[2(m [2nTl . 

2. Comparison to other nearby populations with pub- 
lished Lya measurements (iScarlata et al.l l2009t 
Cowie etldl l20Ill: [Haves et al.l 12014 lOstlin et al.l 
2014 : iMartin et al.f 201511 . sh^ed that the Green 
Peas have (on average) lower masses, higher SFRs 
and specific SFRs, brighter Muv, bluer UV slopes, 
higher Who, lower dust extinction, and smaller 
sizes. 

3. Despite similarly low dust extinction in the Green 

Pea sample, we measure Lya escape fractions 
(within the COS aperture) that span a factor of 
ten: = 0.05 - 0.62. 
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Fig. 15. — The wings and peaks of the Lyo: line profiles are 
compared to the absorption line profiles. The Lya line (blue, with 
peaks not shown), exhibits blueshifted emission which extends to 
the same maximal velocity as the Ly/3 (or Ly 7 ) lines (shown in 
black). 

4. We detect LIS absorption in 8/10 Green Peas, 

and confirm previous findings that weaker Wlis 
is associated with stronger Lya emission (e. g. 
iShaolev et al.l[200^ . albeit, extending prior mea¬ 
surements to higher WLya and weaker Wlis- Al¬ 
though the correlation between Wlis and WLya 
is not significant in our data, a tentative correla¬ 
tion is found when we use directly, instead 

of W Lya as a proxy. Furthermore, the robustness 
of this relation is also supported by detection in 
two independent measures: Si II I260A and C II 
I334A. At the same time, we show that the equiv¬ 
alent width of Si III A1206 does not vary among 
the present sample. 

5. Two of the ten Green Peas show no LIS metal de¬ 

tections, and these are found among those with the 
highest Lya escape fraction, > 0.40, within 

the COS aperture. 

6. Absorption in the H I Lyman series (/3, 7, 5 and/or 
e) is detected in all ten Green Peas. Unlike the 
LIS metal lines, the H I absorption appears uni¬ 
form in kinematics and equivalent width across the 
sample. Furthermore, and in contrast to the LIS 
metal lines, the Lyman series lines are unambigu¬ 


ously saturated, with little to no residual intensity 
at line center. 

7. Kinematic measurements from Si II and C II show 

no correlation between and outflow velocities 
{vc or Vmax) in the LIS metal lines, although two 
of the three highest Green Peas cannot be 

included due to the absence of detectable LIS ab¬ 
sorption. The data do show a marginal correlation 
between and the outflow velocities {vc and 

Vmax) for Si III, but because this line does not trace 
neutral hydrogen, this trend is not compelling evi¬ 
dence that scattering in the cool, neutral phase of 
the outflow is aiding Lya escape. 

8. The Lya escape fraction, shows a signifi¬ 

cant anti-correlation with the velocity separation 
between the blue and red peaks, which is driven by 
shifts in the blue peak velocity. This trend is con¬ 
sistent with a changing neutral hydrogen column 
density, where lower H I density allows more Lya 
to emerge near the systemic velocity. 

9. The velocity of the blue Lya emission spans the 
same range of velocities as the H I absorption lines, 
demonstrating (for the first time with H I) that the 
Lya emission arises from an envelope of gas span¬ 
ning several hundred km s“^. At the same time, the 
maximal outflow velocities reached by the Lyman 
series lines are consistently bluer than the maxi¬ 
mum outflow velocities of the LIS metal lines; the 
good agreement between the velocities probed by 
blue Lya emission and Ly/3/Lyy absorption implies 
that compared to LIS metals, these lines are a bet¬ 
ter proxy for the gas that scatters Lya photons. 

These observations allow us to take a closer look at 
how the interstellar and circumgalactic gas regulates the 
escape of Lya photons from the central few kpc of the 
Green Peas (as probed by the GOS aperture). In past 
studies, a few different Lya escape mechanisms have been 
considered as possibly important. In the remainder of 
this section we take a comprehensive look at the present 
Green Pea data, and explore which Lya escape models 
are consistent with our results. 

Galaxy Outflow Kinematics — Galactic outflows are one 
mechanism for increasing the transmission of Lya pho¬ 
tons. In this case, the Lya photons can scatter in H I 
gas, which is Doppler shifted with respect to the ISM. 
The scattered photons are then out of resonance with the 
ISM and esc ape more easi l y. Pre v ious studies have note d 
this effect: iK unth et all (1199^ . IWofford et al.l (1201311 , 
and IRivera-Thorsen et al.l (1201511 show that Lya emit¬ 
ting galaxies have outflow velocity v ^ —100 to —200 
km s“^, while Lya absorbing galaxies show no outflows 
(LIS absorption is near u ^ 0). On one hand, the Green 
Peas are consistent with these previous measurements; 
all LIS lines are measured with Vc ^ —100 to —200 km 
s“^, and all the Green Peas show Lya emission. On 
the other hand, despite similar outflow velocities, the 
Green Peas show a factor of eight difference in the Lya 
escape fraction (0.05 < < 0.41). Taken together 

with the previous studies, these data suggest that out¬ 
flows may be necessary to permit Lya escape, but they 
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Fig. 16.— A summary of the kin ematic features in the Green Pea data is compared to results from stacked spectra at 2 : ~ 2 — 4 
(ISteidel et al .1120 101 : 1 Jones et al.1120121') . For the Green Peas, multiple transitions are combined to give the range of velocities probed by Ft I 
absorption (green) or LIS metals (Si II and C II; orange and red). In these cases, bars represent the error weighted mean and its uncertainty. 
Additionally, the Green Pea data are sorted in order of increasing (bottom to top). In contrast to the results from high-redshift 

composite spectra (shown at bottom), most of the Green Peas show Lyo emission emerging much closer to the systemic velocity, near Vc 
rather than the LIS Vmax ■ 


are insufficient to explain the wide range of Lya emission 
strengths. Ultimately, more work is needed to clarify the 
role of ontflows, and confirm that the relation between 
outflow velocity and Lya escape is not secondary to some 
other physical characteristic like the ionization state or 
geometry of outflowing gas. 

Neutral Hydrogen Covering — Studies using stacked, com¬ 
posite spectra of high-redshift galaxies (3 < 0 < 4) have 
suggested that Lya escape is regu lated by the cover¬ 
ing fraction in neutral hydrogen, /c (iShaolev et al.ll200l : 
I.Tones et al.ll2012Ll2013[ ). The motivation for this model is 
the correlation between Wlis and WLya- Since the LIS 
metal lines appear saturated in composite spectra, but 
the absorption is not black, the interpretation is that this 
absorbing gas does not cover the entire galaxy. Weaker 
Wl/s suggests lower /c in H I gas, resulting in greater 



While our data confirm the correlation between Wlis 
and W Lya (see (jH) , we cannot attribute this relation to 
non-uniform covering in the H I. By directly probing the 
H I Lyman absorption series, rather than using metals 
as a proxy, we showed that the lines are saturated and 
opaque (or nearly so) in all ten Green Peas (see Figure 
[ini)- Hence, we conclude that low covering of HI cannot 
explain the variations in f^fiN seen in our sample. 

Besides affecting our interpretation of Lya escape. 


the addition of the Lyman series lines requires that we 
re-evaluate our picture of how the CGM is structured 
around high-redshift galaxies. Previous analyses of UV 
spectra support a model where saturated LIS metal ab¬ 
sorption arises from cool, dense clo uds of gas entrained 
in a more highly ion ized outflow (jShaDlev et al.l l2003t 
I Heckman et al.l 1201 iD . For the Green Peas, the LIS 
metal lines are difficult to interpret; because of their low 
S/N , and the probable contribution from emission fill¬ 
ing (|Prochaska et al.ir2flllHScarlata fc Panagial2015f) . we 
cannot determine whether the lines are optically thin or 
optically thick. On one hand, if the LIS metals from 
the Green Peas are optically thick, then the data re¬ 
quire a density dependent covering fraction, fc- This 
scenario requires a pervasive low density, low ionization 
component surrounding higher density clumps; only the 
latter give rise to detectable metal absorption. On the 
other hand, the data could also be described by a ho¬ 
mogeneous low density envelope of gas, which gives rise 
to optically thin LIS metal absorption. The Si II lines 
that we observe become optically thick at line center for 
Nsi II ~0.9 — 3.8x 10^^ cm~^; assumin g solar abun¬ 
dance [Si/H] = -4.5 (jAsplund et al.l[200^ . Z = O.IZq, 
and 100% of the Si exists in Si II, the implied H I columns 
would be around 10^® — 10^® cm^. These are consis¬ 
tent with our observations of the Lyman absorption lines. 
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which we inferred were on the flat part of the curve of 
growth. 

A difference between the spatial distribution of 
detected me t als a nd H I has been noted before. 
IQuider et al.l ()201Clll inferred the existence of a pervasive 
H I component in the lensed “Cosmic Eye” at z ~ 3. 
Similar to the Green Peas, the spectrum of this galaxy 
showed partial covering in metal lines, yet complete cov¬ 
ering in damped Lya absorption. The present COS 
data demonstrate that the “Cosmic Eye” is not alone 
in this characteristic. These observations suggest that a 
widespread, low density H I component may be common. 

Spatially Extended Lya Emission — In Sjd] we noted that 
extended Lya emission, falling outside the COS aperture, 
may be significant for this sample. Three galaxies have 
serendipitous Lya measurements that should approxi¬ 
mate a “total” flux. The LARS imaging showed that 
the COS aperture captured 40% of the flux of 0926-1-4427 
(jHaves et al.ll2014l l , and GALEX grism spectra show that 
we detected 60 and 75% of the flux for 1133-1-6514 and 
1219-1-1526. Intriguingly, for this small subsample, the 
fraction of Lya flux detected by GOS increases with 
However, an aperture correction alone still does 
not account for all of the Lya photons produced in H II 
regions; the small but significant amount of dust mea¬ 
sured in 521 is probably still important, especially when 
coupled with resonant scattering. Systematic uncertain¬ 
ties on the absolute dust correction for Lya preclude a 
full and accurate accounting of the flux. 

The emission and absorption line profiles in our UV 
spectra give another means for testing the significance 
of aperture effects. First, in a qualitative sense, aper¬ 
ture losses can mimic some of the features observed in 
our Lya spectral profiles. For a spherically expanding 
shell, the outermost regions are expanding in the plane 
if the sky, with projected v elocity u ^ 0 (see Figure 2 
in iScarlata fc Panagiall201^ . Consequently, including a 
greater contribution from these regions would increase 
the contribution from Lya emission near the systemic 
velocity, while also increasing At the same time, 

a greater contribution from spatially extended, scattered 
emission would also increase the amount emission from 
metals. This emission would both “fill in” the resonant 
Si and C absorption lines in our data, and would also pro¬ 
duce fluorescent Si II* and C II* emission. In the absence 
of an aperture, and for a spherical geometry, we expect 
the equivalent widths of em ission and absorption t o sum 
to zero: Wats + hEem = 0 (jProchaska et al.l [201 Ih . But 
if the finite COS aperture misses so me of the scattering 
CGM gas, Wem will be decreased (|Scarlata fc Pana^ 
1201511 . We tested whether C II and Si II transitions in 
the present sample were consistent with Wats + Wem = 0, 
but found that this measurement is extremely sensitive 
to the precise placement of the continuum. Because of 
the 10% normalization uncertainty appropriate for our 
data, we cannot determine whether scattered Si II and 
C II emission is missing from the COS aperture. 

Ultimately, direct detection of extended Lya emission, 
similar to the approach adopted for the LARS galaxies 
((Haves et al.l 12013112014 lOstlin et al.l 120141 . is needed. 
Such data would clarify how much total Lya emission 
emerges from the Green Peas, while determining its spa¬ 
tial extent and distribution, and whether any regions 


show Lya absorption. 

Neutral Hydrogen Density — Another galaxy property 
that should impact Lya output is the neutral hydrogen 
column density. With a greater density of H I, the Lya 
photons will undergo increased scattering. In this case, 
they are more susceptible to absorption by dust, and 
must diffuse further in frequency from the line core be¬ 
fore escaping. 

Variations in the H I density offer the most promis¬ 
ing explanation for the range of that we detect in 
the Green Peas. In (17.21 ('see Figure (141). we showed that 
the correlates tightly with the velocity separation 

between the blue and red Lya emission peaks. This cor¬ 
relation is driven mostly by a shift in the blue emission 
peak, which appears closer to the systemic velocity when 
fesc°‘ is high. Since we see no correlation betw een 
and the outflow kinematics discussed in gLU the kine¬ 
matic trends seen in the Lya lines mo re probably arise 
from a sequence in H I density. In fact, iVerhamme et all 
(|2014ll have already attributed small peak separations to 
low H I column densities. In this model, the peaks are 
closer together when low densities allow the Lya pho¬ 
tons to escape close to the systemic velocity. Hence, in 
the present sample, we attribute the increase in to 
a decrease in the neutral hydrogen density. It is plausible 
that decreased dust contents correspond with lower H I 
column density, but more sensitive dust measurements 
are needed to test this effect. Finally, we note that vari¬ 
ations in H I density may explain the correlation between 
feJi°‘ and Wlis which we showed in (JU if the metal ab¬ 
sorption arises from optically thin gas in a homogeneous 
envelope, the lower columns of neutral H I would also 
imply weaker LIS absorption. 

9. CONCLUSIONS 

In this paper, we have presented a COS/FUV spec¬ 
troscopic study of ten Green Pea galaxies at z ~ 0.2. 
We have focused on understanding how Lya photons es¬ 
cape from the ISM and CGM. The Lya emission line 
is frequently used in high-redshift studies, not only to 
discover objects, but also to infer crude constraints on 
the properties of galaxies and the IGM. Yet our limited 
understanding of how Lya photons escape from galaxies 
hampers our ability to fully interpret high-redshift obser¬ 
vations (|Diikstra et al.l[Ml4 iDressler et al.l 120151 1. The 
detection of Lya emisison in all ten of the Green Peas 
suggests that in the lowest mass, least dusty, highest SFR 
galaxies this feature may appear in the majority of cases. 

Beyond simple detection of Lya, the UV absorption 
line data provides important constraints on the ISM and 
CGM gas that regulate Lya output. These data show 
that while cool outflows traced by LIS metal absorption 
may play a role in permitting Lya escape, they do not 
explain the widely varying Lya strengths observed in 
the Green Peas. Moreover, aided by the low redshifts 
of these galaxies and the absence of Lya forest absorp¬ 
tion, for the first time we are able to directly analyze 
H I absorption in the Lyman series. These transitions 
show clear differences from the LIS metal lines which are 
usually used as a proxy for H I; most importantly, the 
H I covering fraction, /c, is uniformly high for gas with 
Nh^ > 10^® cm“^. The Lya photons do not escape the 
Green Peas through holes which are completely devoid of 
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CGM/ISM gas. Instead, the kinematic variations in the 
Lya emission profiles are strongly suggestive that H I col¬ 
umn density is the primary characteristic regulating Lya 
escape in these galaxies. The ubiquity of double-peaked 
prohles, and the small velocity separation between the 
peaks strongly suggests that a low column density of H I 
gas is a typical characteristic among Green Pea galaxies. 

This study demonstrates that nearby, high-redshift 
analog galaxies are useful local laboratories, where a 
wealth of high fidelity data can provide challenges for 
frequently adopted models. By observing the “Green 
Pea” galaxies, we have provided a fresh look at objects 
which, to our knowledge, are excellent analogs for young 
galaxies in an early Universe. Yet much work remains 
before we can be certain that the conclusions drawn here 
are applicable at moderate to high-redshifts. Gertainly, 
spectroscopy with James Webb Space Telescope will clar¬ 
ify the overlap between high-z galaxies and their nearby 
counterparts. Ultimately, this study serves as a strong 
motivation for future investigations. An even greater 


understanding of these galaxies could be achieved with 
maps of the spatially extended Lya emission, as well as 
sensitive measures of the small amount of dust present 
in these systems. 
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APPENDIX 

COS SPECTRA OF THE GREEN PEAS 
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Fig. A1.— COS G130M + G160M spectrum of the Green Pea galaxy 0303-0759. Gaps in the spectrum demarcate the wavelengths 
between the COS FUV A and B segments. Red markers indicate the expected locations of ISM features in the Green Pea and blue labels 
show stellar features. Black vertical lines indicate Milky Way features, whereas 0 symbols mark geocoronal emission. This sight line 
contains damped Milky Way Lyo; absorption, apparent around Xrest 1042 A. 
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Fig. A2.— Same as Figure IaTI but for 1244+0216. This sight line does not contain damped Milky Way Lya absorption. 



Fig. A3.— Same as Figure [All but for 1054+5238. 


Fig. A4.— 



Same as Figure \Kl\ but for 1137+3524. 
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Fig. A7.— Same as Figure iah but for 1424+4427 . This spectrum has shorter wavelength coverage because the G160M observation 
failed and our program did not qualify for a repeat. 
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